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Abstract 

We derive predictions for the Dirac phase 8 present in the 3x3 unitary neutrino mixing 
matrix U = U^U^j, where [4 and U^, are 3x3 unitary matrices which arise from the 
diagonalisation, respectively, of the charged lepton and the neutrino mass matrices. We 
consider forms of Ue and U,, allowing us to express <5 as a function of three neutrino 
mixing angles, present in U, and the angles contained in U,,. We consider several forms 
of Ui, determined by, or associated with, symmetries, tri-bimaximal, bimaximal, etc., for 
which the angles in [7^ are fixed. For each of these forms and forms of allowing one to 
reproduce the measured values of the neutrino mixing angles, we construct the likelihood 
function for cos 6, using i) the latest results of the global fit analysis of neutrino oscillation 
data, and ii) the prospective sensitivities on the neutrino mixing angles. Our results, in 
particular, confirm the conclusion, reached in earlier similar studies, that the measurement 
of the Dirac phase in the neutrino mixing matrix, together with an improvement of the 
precision on the mixing angles, can provide unique information as regards the possible 
existence of symmetry in the lepton sector. 
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1 Introduction 


Understanding the origin of the observed pattern of neutrino mixing, establishing the status 
of the CP symmetry in the lepton sector, determining the type of spectrum the neutrino 
masses obey and determining the nature — Dirac or Majorana — of massive neutrinos are 
among the highest priority goals of the programme of future research in neutrino physics (see, 
e.g., [^). One of the major experimental efforts within this programme will be dedicated to 
the searches for CP-violating effects in neutrino oscillations (see, e.g., [^[^). In the reference 
three neutrino mixing scheme with three light massive neutrinos we are going to consider 
(see, e.g., [^), the CP-violating effects in neutrino oscillations can be caused, as is well 
known, by the Dirac CP violation (CPV) phase present in the Pontecorvo, Maki, Nakagawa, 
Sakata (PMNS) neutrino mixing matrix. Predictions for the Dirac CPV phase in the lepton 
sector can be, and were, obtained, in particular, combining the phenomenological approach, 
developed in [4j|^ and further exploited in various versions by many authors with the aim 
of understanding the pattern of neutrino mixing emerging from the data (see, e.g., [^13|), 
with symmetry considerations. In this approach one exploits the fact that the PMNS mixing 
matrix U has the form [^: 

U = UlU, = {Ue)^^UuQo: (1) 

where Ue and U^, are 3x3 unitary matrices originating from the diagonalisation, respectively, 
of the charged lepton and neutrino mass matrices. In eq. Q Ue and are CKM-like 3x3 
unitary matrices, and 'k and Qq are diagonal phase matrices each containing in the general 
case two physical CPV phases 0 

^'= diag , Qo = diag(l,U^ . (2) 

It is further assumed that, up to subleading perturbative corrections (and phase matrices), 
the PMNS matrix U has a specific known form Ui^ that is dictated by continuous and/or 
discrete symmetries, or by arguments related to symmetries. This assumption seems very 
natural in view of the observation that the measured values of the three neutrino mixing 
angles differ from certain possible symmetry values by subdominant corrections. Indeed, the 
best fit values and the 3 (t allowed ranges of the three neutrino mixing parameters sin^ 0 i 2 , 
sin^ 023 and sin^ 0i3 in the standard parametrisation of the PMNS matrix (see, e.g., [^), 
derived in the global analysis of the neutrino oscillation data performed in read 

(sin^ 0 i2)bf = 0.308 , 0.259 < sin^ 9^ < 0.359 , (3) 

(sin^ 023 )bf = 0.437 (0.455), 0.374 (0.380) < sin^ 023 < 0.626 (0.641), (4) 

(sin^ 0 i3)bf = 0.0234 (0.0240), 0.0176 (0.0178) < sin^ 0i3 < 0.0295 (0.0298), (5) 

where the value (the value in parentheses) corresponds to > 0 < 0), i.e., 

neutrino mass spectrum with normal (inverted) ordering]^ (see, e.g., [^). In terms of angles, 
the best fit values quoted above imply; 0i2 = 7r/5.34, 0i3 = 7r/20 and 023 = 7r/4.35. Thus, for 

^If the charged lepton mass term is written in the right-left convention, the matrix Ue diagonalises the 
hermitian matrix mI^Me, UIM^^MeUb = diag(rrie, rn,/, m^), Me being the charged lepton ma ss matrix. 

^The phases in the matrix Qo contribute to the Majorana phases in the PMNS matrix lid]. 

^Similar results were obtained in the global analysis of the neutrino oscillation data performed in nii|. 
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instance, O 12 deviates from the possible symmetry value vr/4, corresponding to the bimaximal 
17,^, by approximately 0.2, ^13 deviates from 0 (or from 0.32) by approximately 


mixing 

0.16 and 023 deviates from the symmetry value vr/d by approximately 0.06, where we used 
sin^ 023 = 0.437. 

Widely discussed symmetry forms of tJu include: i) tri-bimaximal (TBM) form js 
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bimaximal (BM) form, or due to a symmetry corresponding to the conservation of the lepton 
charge L' = L^ — L^ — (LC) [iTp^, hi) golden ratio type A (GRA) form 20,^, iv) golden 


ratio type B (GRB) form [^, and v) hexagonal (HG) form [13[[^ . For all these forms the 
matrix represents a product of two orthogonal matrices describing rotations in the 1-2 and 
2-3 planes on fixed angles 


^2 and 02^3: 


U, = R2z{e'i^) Ri2{e\2), 


where 


1 

( COS 05^2 

sin 0^2 


1 1 


0 

0 \ 

Ri 2 (0/2) = 

— sin 05^2 

cos 05(2 

0 

, 7?23 ( 02 " 3 ) = 

0 

cos 023 

sin 0^3 

\ 

V 0 

0 


' \ 

Vo 

-sin 0^3 

cos 0^3/ 


( 6 ) 


(7) 


Thus, Ui, does not include a rotation in the 1-3 plane, i.e., 61 ^ = 0. Moreover, for all the 
symmetry forms quoted above one has also 023 = — 7r/4. The forms differ by the value of 
the angle 0)^2) and, correspondingly, of sin^ 0 i 2 '- for the TBM, BM (LG), GRA, GRB and HG 
forms we have, respectively, sin^ 05^2 = 1/3, 1/2, (2 -|- r)~^ = 0.276, (3 — r)/4 = 0.345, and 
1/4, r being the golden ratio, r = (1 -|- \/5)/2. 

As is clear from the preceding discussion, the values of the angles in the matrix Un, which 
are fixed by symmetry arguments, typically differ from the values determined experimentally 
by relatively small perturbative corrections. In the approach we are following, the requisite 
corrections are provided by the angles in the matrix Ue- The matrix Ue in the general case 
depends on three angles and one phase [^. However, in a class of theories of (lepton) flavour 
and neutrino mass generation, based on a GUT and/or a discrete symmetry (see, e.g., [24[|^ ), 
Ue is an orthogonal matrix which describes one rotation in the 1-2 plane, 


Ue = Ri^{ei2 ), 


( 8 ) 


or two rotations in the planes 1-2 and 2-3, 

Ue = R2^{9l^)R^i{el2), 


(9) 


0^2 s-iid 0^3 being the corresponding rotation angles. Other possibilities include Ue being an 
orthogonal matrix which describes i) one rotation in the 1-3 plane 

Ue = R-,^{e%). (10) 

or ii) two rotations in any other two of the three planes, e.g., 

Ue = R 2 ^{ei^)R^^{ei^), or (11) 

Ue = R-^^{6l^)R^i{ei2). ( 12 ) 

■^The case of representing a rotation in the 2-3 plane is ruled ont for the five symmetry forms of W listed 
above, since in this case a realistic value of @13 7 ^ 0 cannot be generated. 
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( 12 ) is a matter of convenience — this allows us 


The use of the inverse matrices in eqs. 
to lighten the notations in expressions which will appear further in the text. 
It was shown in (see also 


31 


) that for Up and Ue given in eqs. ([^ and ([9|), the Dirac 
phase 5 present in the PMNS matrix satisfies a sum rule by which it is expressed in terms 
of the three neutrino mixing angles measured in the neutrino oscillation experiments and the 
angle 0 i2- In the standard parametrisation of the PMNS matrix (see, e.g., [^) the sum rule 
reads 


30 


cos (5 = 


tan 023 


sin 2012 sin 013 


[cos 2012 + Gi 2 — cos^ 0 i 2 ) (l — cot^ 023 sin^ 013 )] . (13) 


For the specific values of 0i2 = vr/d and 0i2 = sin ^(l/v/S), i.e., for the BM (LC) and TBM 
forms of Up, eq. (13) reduces to the expressions for cos J derived first in 
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On the basis of 

the analysis performed and the results obtained using the best fit values of sin^ 0 i 2 , sin^ 0 i 3 


and sin^ 023 , it was concluded i 


m 
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in particular, that the measurement of cos <5 can allow 


one to distinguish between the different symmetry forms of the matrix Up considered. 
Within the approach employed, the expression for cos <5 given in eq. (13) is exact. In 


the correction to the sum rule eq. (13) due to a non-zero angle 0f3 <C 1 in Ue, corresponding 
to 

s)^r3(0?3)^r2'(012) (14) 


Ue = i?23 (0|3i 


with I sin 0 ^ 3 1 <C 1, was also derived. 

Using the best fit values of the neutrino mixing parameters sin^ 0i2, sin^ 0i3 and sin^ 023, 
found in the global analysis in predictions for cos J, 6 and the rephasing invariant 


1 


Jcp =lm{UeiU* 3 Ue 3 U^i} = - sinsin 2013 sin 2023 sin2012 cos 013 , 


(15) 


which controls the magnitude of CP-violating effects in neutrino oscillations [^, were pre¬ 
fer each of the five symmetry forms of Up — TBM, BM (LC), GRA, GRB and 
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sented in 


HG — considered. 


Statistical analysis of the sum rule eq. (13) predictions for 6 and Jcp (for cos 5) using 
the current (the prospective) uncertainties in the determination of the three neutrino mixing 
parameters, sin^ 0 i 2 , sin^ 0 i 3 , sin^ 023 , and S (sin^ 0 i 2 , sin^ 0 i 3 and sin^ 023 ), was performed 
for the five symmetry forms — BM (LC), TBM, GRA, GRB and HG — of Up. Using 


m 
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the current uncertainties in the measured values of sin^ 0 i 2 , sin^ 0 i 3 , sin^ 023 and it was 
found, in particular, that for the TBM, GRA, GRB and HG forms, Jcp 7 ^ 0 at 5 cj, 4cr, 4 ct 
and 3cj, respectively. For all these four forms | Jcp| is predicted at Sa to lie in the following 
narrow interval 34 : 0.020 < | Jcp| < 0.039. As a consequence, in all these cases the CP- 


violating effects in neutrino oscillations are predicted to be relatively large. In contrast, for 
the BM (LC) form, the predicted best fit value is Jcp = 0, and the CP-violating effects 
in neutrino oscillations can be strongly suppressed. The statistical analysis of the sum rule 
predictions for cos (5, performed in 34 by employing prospective uncertainties of 0.7%, 3% 


and 5% in the determination of sim0i2, sin^ 0i3 and sin^ 023, revealed that with precision 
in the measurement of 6 , A6 = (12° — 16°), which is planned to be achieved in the future 
neutrino experiments like T2HK and ESSi^SB [^, it will be possible to distinguish at 3a 

®We would like to note that the recent statistical analyses performed in [l 5 [|l 6 | showed indications/hints 
that 5 = 37r/2. As for sin^ 612, sin^ 613 and sin^ 623, in the case of 5 we utilise as “data” the results obtained 
in ref. 
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between the BM (LC), TBM/GRB and GRA/HG forms of tJu- Distinguishing between the 
TBM and GRB forms, and between the GRA and HG forms, requires a measurement of 6 
with an uncertainty of a few degrees. 

In the present article we derive new sum rules for cos <5 using the general approach em¬ 
ployed, in particular, in [^[^. We perform a systematic study of the forms of the matrices 
Ue and Uu, for which it is possible to derive sum rules for cos <5 of the type of eq. (13), but for 


which the sum rules of interest do not exist in the literature. More specifically, we consider 
the following forms of Ue and U^: 

A. Uy = i?23(023)-^i2(^i2) with 023 = — 7 r /4 and 6'(2 corresponding to the TBM, BM (LG), 


GRA, GRB and HG mixing, and i) 


Ue = R^^{e' 
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ii) Ue — R2Z (^23 


)R 


13^ (^13) 


and hi) 


Ue = R^i{ei^)Rl^{ei2)-, 


B. lly = i?23(^23)-^i3(^i3)'^i2(^i2) with 023; ^13 ^nd 0^12 fixed by arguments associated with 
symmetries, and iv) Ue = Ri 2 {'Si 2 )i s-'^d v) Ue = R^ 3 ^( 0 f 3 ). 

In each of these cases we obtain the respective sum rule for cosd. This is done first for 
023 = — 7 r /4 in the cases listed in point A, and for the specific values of (some of) the angles 
in Uy^ characterising the cases listed in point B. For each of the cases listed in points A and B 
we derive also generalised sum rules for cos 5 for arbitrary fixed values of all angles contained 
in Uu (i.e., without setting 023 = — 7 r /4 in the cases listed in point A, etc.). Next we derive 
predictions for cos (5 and Jcp (cos (5), performing a statistical analysis using the current (the 
prospective) uncertainties in the determination of the neutrino mixing parameters sin^ 0 i 2 , 
sin^ 013 , sin^ 023 and 6 (sin^ 0 i 2 , sin^ 0 i 3 and sin^ 023 ). 

It should be noted that the approach to understanding the experimentally determined 
pattern of lepton mixing and to obtaining predictions for cos <5 and Jcp employed in the 
present work and in the earlier related studies 
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and , is by no means unique — it is one 
of a number of approaches discussed in the literature on the problem (see, e.g., [35}|37] ). It is 
used in a large number of phenomenological studies (see, e.g., muioHilPi ) as well as in 
a class of models (see 


24 -29 ,^) of neutrino mixing based on discrete symmetries. However, 


it should be clear that the conditions of the validity of the approach employed in the present 
work are not fulfilled in all theories with discrete flavour symmetries. For example, they are 
not fulfilled in the theories with discrete flavour symmetry A( 6 n^ 
the S '4 flavour symmetry constructed in 
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studied in 

and in the models discussed in (4^ . 


40 
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with 
Further, the 

conditions of our analysis are also not fulfilled in the phenomenological approach developed 
and exploited in 36,37 . In these articles, in particular, the matrices Ue and Uy are assumed 


to have specific given fixed forms, in which all three mixing angles in each of the two matrices 
are fixed to some numerical values, typically, but not only, 7 r/ 4 , or some integer powers n of 
the parameter e = Oc, Oc being the Gabibbo angle. The angles = (Oc)^ with n > 2 are 
set to zero. For example, in 37 the following sets of values of the angles in Ue and Uy have 


been used: ( 0 i 2 , 6 'f 3 > ^ 23 > ^ 13 > ^ 23 ) = (*, vr/4, 7 r/ 4 , =t=, 7 r/ 4 , =t=) and (*, =t=, 7 r/ 4 , 7 r/ 4 , *, *), where 

means angles not exceeding 9c- None of these sets correspond to the cases studied by us. 
As a consequence, the sum rules for cos (5 derived in our work and in are very different. 
In 37 the authors have also considered specific textures of the neutrino Majorana mass 
matrix leading to the two sets of values of the angles in Ue and Uy quoted above. However, 
these textures lead to values of sin^ 023 or of sin^ 0 i 2 which are strongly disfavoured by the 
current data. Although in a large variety of forms of Ue and Uy have been investigated, 
none of them corresponds to the forms studied by us, as can be inferred from the results on 
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the values of the PMNS angles 0 i 2 , ^13 and @23 obtained in and summarised in Table 2 
in each of the two articles we have cited in 36 . 

Our article is organised as follows. In Section we consider the models which contain 
one rotation from the charged lepton sector, i.e., IJe = = -^r 3 ^(^f 3 )) and two 

rotations from the neutrino sector: Uu = .^ 23 (^ 23 ) ^ 12 (^ 12 )- In these cases the PMNS matrix 
reads 


u = RiMj) ^ ^ 23 (^ 23 ) Ri2{e'i2) Qo , ( 16 ) 

with (ij) = (12), (13). The matrix Ui, is assumed to have the following symmetry forms: TBM, 
BM (LC), GRA, GRB and HG. As we have already noted, for all these forms 023 = 
but we discuss also the general case of an arbitrary fixed value of 023- The forms listed above 
differ by the value of the angle 0 ^ 2 ) which for each of the forms of interest was given earlier. 
In Section we analyse the models which contain two rotations from the charged lepton 
sector, i.e., Ue = ^^ 3 ^( 023 ) -I^r 3 ^(^i 3 )! or Re = -Rr 3 ^(^i 3 ) I^r 2 ^(^i 2 )) and|^two rotations from the 
neutrino sector, i.e.. 


U = R,,{et^) RkMi) T ^23(02^3) ^ 12 ( 0 ^ 2 ) Qo , (17) 

with (ij) — {kl) = (13) — (23), (12) — (13). First we assume the angle 023 to correspond to the 
TBM, BM (LC), GRA, GRB and HG symmetry forms of Uu- After that we give the formulae 
for an arbitrary fixed value of this angle. Further, in Section we generalise the schemes 
considered in Sectionby allowing also a third rotation matrix to be present in ily-. 


U = R^j) 


T R23( 02 ^ 3 ) 7^13(0^3) 7?12 (05(2) Qo 


(18) 


with {ij) = (12), (13), (23). 

Using the sum rules for cos 5 derived in Sections in Section we obtain predictions 

for cos d, 5 and Jcp for each of the models considered in the preceding sections. Section 
contains summary of the results of the present study and conclusions. 

We note finally that the titles of Sections 2-4 and of their subsections reflect the rotations 
contained in the corresponding parametrisation, eqs. (16) - (18). 


2 The Cases of 9R — 612 ) Rotations 


In this section we derive the sum rules for cosd of interest in the case when the matrix 
Uu = 7?23(023) 77i2(0i 2) with fixed (e.g., symmetry) values of the angles 023 and 05 ^ 2 , gets 
correction only due to one rotation from the charged lepton sector. The neutrino mixing 
matrix U has the form given in eq. (|16[). We do not consider the cases of eq. (16) i) with 


{ij) = (23), because the reactor angle 0i3 does not get corrected and remains zero, and ii) 


with {ij) = (12) and 023 = —tt/A, which has been already analysed in detail in 30,34 . 


2.1 


The Scheme with 0®2 


(023, 012) Rotations 


For 0-3 = —7r/4 the sum rule for cos <5 in this case was derived in ref. and is given in 


eq. (50) therein. Here we consider the case of an arbitrary fixed value of the angle 0 ^ 3 . Using 


31 


° We consider only the “standard” ordering of t he two rotations in Ue, see 
i^ 23 {^ 23 ) Rr 2 ^(^i 2 ) has been analysed in detail in 30 31 34 and will not be discussed by us. 


The case with Ue = 
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eq. (16) with (ij) = ( 12 ), one finds the following expressions for the mixing angles @13 and 
023 of the standard parametrisation of the PMNS matrix: 


sin^ 013 = |17e3p = sin^ 0i2 sin^ 0. 


sin^ 023 = 


\U, 


/iSl 


23 ) 

sin^ 023 — sin^ 0 


13 


1 - \Ue3\^ 


1 — sin^ 013 


(19) 

( 20 ) 


Althongh eq. (13) was derived in 30 for 0^3 = —vr/d and Ue = i?23^(0|3)i?)”2^(0f2), it is not 
difficult to convince oneself that it holds also in the case under discussion for an arbitrary 
fixed value of 023- The sum rule for cos <5 of interest, expressed in terms of the angles 0 i 2 , 
013, 0^2 and 023, can be obtained from eq. ( |l3[ ) by using the expression for sin^ 623 given in 
eq. (20). The result reads: 


cos (5 = 


(cos 2013 ~ cos 2023 ) ^ 
\/2 sin 2012 sin 0 i 31 cos 0 : 


cos 20 


12 


+ (sin^ 012 — cos^ 012 ) 


231 L 

.2 


2 sin^ 023 — (3 + cos 20 ^ 3 ) sin^ 0 i 3 


cos 2013 ~ cos 20 ; 


23 


( 21 ) 


Setting 023 = — 7 r /4 in (21), one reproduces the sum rule given in eq. (50) in ref. 30 


2.2 The Scheme with 0^3 — (023,012) Rotations 

In the present subsection we consider the parametrisation of the neutrino mixing matrix given 


in eq. (16) with (ij) = (13). In this set-up the phase 1 /^ in the matrix ih is unphysical (it can 


be absorbed in the field) and therefore effectively T = diag(l,l,e *‘^). Using eq. (16) 
with (ij) = (13) and 623 = —vr/d and the standard parametrisation of U, we get 


sin^ 013 = |Ue3|^ = ^ sin^ 013 , 


sim 023 = 




sin^ 012 = 


l-|C/e3|2 2(1-sin2 0i3) ’ 

|C/e2p 1 


1 - |t^e3p 1 - sin^ 013 [2 


- sin^ 033 cos^ 0v 


^12 


1 


-I- cos 013 sin 032 H— 1 = sin 20^3 cos io sin 032 cos 032 

v2 


( 22 ) 

(23) 

(24) 


From eqs. (22) and (24) we obtain an expression for cos a; in terms of the measured mixing 

sin^ 013 


angles 0 i 2 , 0 i 3 and the known 032 : 

1 — sin^ 033 


COSCJ = 


sin 2032 sin 033(1 — 2 sin^ 033)2 - 


sin^ 032 — sin^ 032 — cos 2032 


1 — sin^ 033 


(25) 


Further, one can find a relation between sin <5 (cos <5) and sincj (cosw) by comparing the 
imaginary (the real) part of the quantity U*iU*^Ue 3 U^i, written by using eq. (16) with (ij) = 


'We note that the expression for cos uj we have obtained coincides with that for cos (j) in the set-up with the 


(ij) = (12) rotation in the charged lepton sector (cf. eq. (46) in 30 ). 
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(13) and in the standard parametrisation of U. For the relation between sin (5 and sinw we 
get 

sin 20^0 , . 

12 , . (26) 


sin 5 = — - 


■ sin ui. 


sin 2012 

The sum rule for cos <5 of interest can be derived by substituting cosw from eq. ( |25[ ) in the 
relation between cos 5 and coscj (which is not difficult to derive and we do not give). We 
obtain 


cos (5 = — 


(1 — 2sin^ 013)2 
sin 2012 sin 0 i 3 


cos 2012 + ^12 — cos^ 0 i 2 ) 


1 — 3 sin^ 013 
1 — 2 sin^ 013 


(27) 


We note that the expression for cos 5 thus found differs only by an overall minus sign from 
the analogous expression for cos 5 derived in |30| in the case of (ij) = (12) rotation in the 
charged lepton sector (see eq. (50) in (30]). 


In eq. (15) we have given the expression for the rephasing invariant Jcp in the standard 
parametrisation of the PMNS matrix. Below and in the next sections we give for completeness 
also the expressions of the Jcp factor in terms of the independent parameters of the set-up 
considered. In terms of the parameters oo, and 0i2 of the set-up discussed in the present 
subsection, Jcp is given by 


Jcp = — 


1 


8^2 


sin to sin 20 i 3 sin 20 i 2 . 


(28) 


In the case of an arbitrary fixed value of the angle 023 the expressions for the mixing 
angles 0 i 3 and 023 take the form 


sin^ 013 = I Ue3 1 ^ = sin^ cos^ 023 , 
\U^3\^ 

1 - |C/e3p 


sm ^23 — 


sin^ 0^3 


1 — sin"^ 013 


(29) 

(30) 


The sum rule for cos <5 in this case can be obtained with a simpler procedure, namely, by 
using the expressions for the absolute value of the element I/^i of the PMNS matrix in the 
two parametrisations employed in the present subsection: 


I U^i I = I cos 023 sin 012 -|- cos 0i2 sin 0i3 sin 623 \ = | cos 023 sin 05^2 1 , 
From eq. © we get 

(cos 2013 + cos 2023 ) ^ 


(31) 


cos 5 = — 


\f2 sin 2012 sin 0 i 3 1 sin 0: 


-|- (sin^ 012 — cos^ 012 ) 


V 

231 L 

Z.2 Qv 


COS 2012 


2 COS^ 023 — (3 — COS 2023) ^13 

COS 2013 + COS 2023 


(32) 


We will use the sum rules for cos <5 derived in the present and the next two Sections to 
obtain predictions for cos <5, 8 and for the Jcp factor in Section 
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3 The Cases of {O^j^Oh) — (6*23,6*12) Rotations 

As we have seen in the preceding Section, in the case of one rotation from the charged lepton 
sector and for ^23 = the mixing angle 623 cannot deviate significantly from 7r/4 due to 

the smallness of the angle 613. If the matrix Uu has one of the symmetry forms considered 
in this study, the matrix Ue has to contain at least two rotations in order to be possible to 
reproduce the current best fit values of the neutrino mixing parameters, quoted in eqs. ^ ~ 
([^ . This conclusion will remain valid if higher precision measurements of sin^ 623 confirm 
that 023 deviates significantly from 7 r/ 4 . In what follows we investigate different combinations 
of two rotations from the charged lepton sector and derive a sum rule for cos <5 in each set-up. 
We will not consider the case (012,^23) “ (^23, ^12), because it has been thoroughly analysed 
in refs. [ 30 , 31,34 , and, as we have already noted, the resulting sum rule eq. ( 13 ) derived 


in 30 holds for an arbitrary fixed value of 62^. 


3.1 The Scheme with ( 01 ^, 91 ^) — (^23, ^12) Rotations 

Following the method used in ref. [^, the PMNS matrix U from eq. © with (ij) 
( 13 ) — ( 23 ), can be cast in the form: 

U = Rn{ 9 t 3 ) Pi R 23023 ) Rl2{e'i2) Q , 
where the angle 623 is determined i) for 9 ^^ = — 7 r /4 by 

sin^ §23 = - (1 — sin 202^ cos(a; — V')) , 


and ii) for an arbitrary fixed value of @23 by 


1 


sin"^ 023 = sin"^ $23 cos^ 023 + cos"^ 0^3 sin"^ 023 + - sin 20^3 sin 2023 cos(cu — V’) • 

The phase matrices Pi and Q have the form: 

Pi = diag(l, 1, e“*“), and Q = QiQo, with Qi = diag(l, 1, e*^), 
where the phases a and /? are given by 

a = 7 + ^/> + w , with 7 = arg (cos 0^3 sin 023 + sin 0^3 cos 0, 


V 

23 ] 5 


/? = 7 — (/), where (j) = arg cos 0^3 cos 9^3 — e sin 0^3 sin 023^ 


Using eq. ( 33 ) and the standard parametrisation of U, we find: 
sin^ 013 = |Pe3p = sin^ ^*13 cos^ §23 , 


irr |2 „;„2 

\U, 


sm U23 = 


sm ^12 = 


>31 


sin"' 023 


1 - 1 - sin^ 013 ’ 

|Pe2p 1 


1 - |Pe3p 1 - sin^ 013 


COS^ 033 sin^ 032 


— - sin 023 sin 2033 sin 2032 cos a + cos^ 032 sin^ 033 sin^ §23 


ikl) = 

( 33 ) 

( 34 ) 

( 35 ) 

( 36 ) 

( 37 ) 

( 38 ) 

( 39 ) 

( 40 ) 


( 41 ) 
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The first two equations allow one to express and 023 in terms of 0i3 and 023- Eq- (41) 
allows us to find cos a as a function of the PMNS mixing angles 0i2, 0i3, 023 and the angle 
0l2- 

sin^ 05^2 cos^ 023 + cos^ 6\2 sin^ ^23 sin^ 0i3 — sin^ 0i2 (l — sin^ 023 cos^ 013 ) 

- ^iTi90^ ^m90 ^m0 -• 

sm ^ 0^2 sin ^c 723 sin (713 

The relation [^between sin 6 (cos 6) and sin a (cos a) can be found by comparing the imaginary 
(the real) part of the quantity U*iU*^Ue 3 U^i, written using eq. (33) and using the standard 
parametrisation of U: 


. sm 201^2 

sm 6 = -^ 

sm20i2 


sina, 


(43) 


^ sm20i9 

cos 5 = -^ 

sin 2012 


U g* 0 

cos a -^tan 622, (cos 20i2 + cos 205^2) 


(44) 


sin 2012 

The sum rule expression for cos (5 as a function of the mixing angles 0i2, 0i3, 023 and 9 12 -, with 
05^2 having an arbitrary fixed value, reads: 

cot 023 


cos 5 = — 


sin 2012 sin 013 


[cos 205^2 + (sin^ 012 — cos^ 012 ) (l — tan^ 023 sin^ 013 )] 


(45) 


This sum rule for cos <5 can be obtained formally from the r.h.s. of eq. (13) by interchanging 
tan023 and cot023 and by multiplying it by (—1). Thus, in the case of 023 = 7 r/ 4 , the 
predictions for cos <5 in the case under consideration will differ from those obtained using 
eq. (13) only by a sign. We would like to emphasise that, as the sum rule in eq. (13), the sum 
rule in eq. (45) is valid for any fixed value of 0 ^ 3 . 

The Jcp factor has the following form in the parametrisation of the PMNS matrix em¬ 
ployed in the present subsection: 


Jcp = - sin 20^3 sin 29^2 sin 2023 cos 023 sin a . 


(46) 


3.2 The Scheme with (6'^25^i3) ~ (^231^12) Rotations 

In this subsection we consider the parametrisation of the matrix U defined in eq. (0 with 
[ij) — (kl) = (12) — (13) under the assumption of vanishing a;, i.e., T = diag(l, 1). In the 
case of non-fixed oj it is impossible to express cos 5 only in terms of the independent angles 
of the scheme. We will comment more on this case later. 

Using the parametrisation given in eq. 0 with 023 = —7r/4 and a; = 0 and the standard 
one, we hnd: 


sm ^13 


sin^ 023 = .. 1 2 = —\ X cos^ 9^2 + ^ sin^ 0^2 sin^ 0i: 

l-|(7e3r COS"=013[2 2 


•2/3 r^e2 


Ue 2 p _ c sin^ 0^2 + g 


1 — sim 013 


(47) 

(48) 

(49) 


30 


“We note that the expression ( |42[ ) for cos a can be obtained formally from the r.h.s. of the eq. (22) for 
by substituting sin 623 with cos 623 and vice versa and by changing its overall sign. 
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where 


= 2 ^^12 S™ ^13 CO® V' ) 

C = cos^ cos 26 i 2 H-^ sin 29'(2 cot ^13 (3 cos 2^33 — 1) 


4 V 2 


^ = cos"^ @33 sin^ 012 + 2 (cO® 2013 — cos 20^3) cos"^ 012 

H-^ sin 205^2 (3 cos 01 ^ sin Of^ — 2 cot 0 ^^ sin^ @13). 

2y2 


( 50 ) 

( 51 ) 

( 52 ) 


The dependence on cos V' in eq. ( 49 ) has been eliminated by solving eq. ( 47 ) for X^. It follows 
from eqs. ( 47 ) and ( 48 ) that sin^ 0 f^ is a function of the known mixing angles @13 and 023 - 


Sin^ 033 = 1 — 2 COS^ 013 COS^ 023 . 

Inverting the formula for sin^ 0i2 allows us to find sin^ 0f2, which is given by 
sin^ 012 = 4 [cos 2012(cos ^^13 + sin^ 013) — cos 20i2 cos^ 013] tan 0i3 + V 2 sin 20 i. 


(53) 


12 


X (3 cos 2013 “ 2 cos 2013 — 1) 


4 cos 2012 ®™ 2013 + V^(3 cos 20i3 — 1) sin 20 


1 -1 


( 54 ) 


Using eqs. ( 47 ) and ( 54 ) we can write cos ip in terms of the standard parametrisation mixing 
angles and the known 0i3 and 0i2: 

sin^ 012 + cos^ 0f2 sin^ 0i3 — 2 sin^ 0i3 


cos tp = 


sin 2012 ®i^^ ^13 


( 55 ) 


We find the relation between sin 5 and sin ip by employing again the standard procedure of 
comparing the expressions of the Jcp factor, Jqp = Im(t/git/* 3 t/e 3 t/i(i), in the two parametri- 
sations — the standard one and that defined in eq. ( flT] ) (with 023 = —7r/4 and oj = 0 ): 

sin 2012 sin ip 


sin (5 = 


4 sin 2012 sin 20i3 sin 023 


2V2 sin 20^3 cos 205(2 + (3 cos 20^=3 - 1) sin 20 


12 


( 56 ) 


where sin20f2 (sin20i3 and 00320^3) can be expressed in terms of 0i2, 0i3, 023 and 0i2 (^13 
and 023) using eq. (Q (eq. ([^). 

We use a much simpler procedure to find cos(i. Namely, we compare the expressions for 
the absolute value of the element Uri of the PMNS matrix in the standard parametrisation 
and in the symmetry related one, eq. with 023 = — 7r/4 and w = 0 , considered in the 
present subsection: 

\Uri\ = I sin 023 sin012 — sin0i3 cos 023 cos0i2e*'^| = | sin0i3 cos0i2 H—^ cos0i3 sin0i2l • ( 57 ) 


^/2 


From the above equation we get for cos 5 : 

. 2 
cos 0 = — 


sin 2012 sin 2023 sin 0i3 


cos^ 023 sin^ 012 sin^ 0i3 + cos^ 0i2 sin^ 023 


— COS^ 013 COS^ 023 cos 0i2 — K \/1 — 2 COS^ 013 COS^ 023 sin 0i2^ 


2 1 


( 58 ) 
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where k = 1 if 6\^ belongs to the first or third quadrant, and k = —1 if is in the second or 
the fourth one. In the parametrisation under discussion, eq. with (ij) — (kl) = (12) — (13), 
^23 = —vr/4 and w = 0, we have: 


Jcp = ^ cos sin 20^2 (^2-\/2 cos 205^2 sin 20^3 + (3 cos 29^^ — 1) sin 20^2^ sin . (59) 


In the case of non-vanishing uj, using the same method and eq. (53), which also holds for 
w 7 ^ 0 , allows us to show that cos 5 is a function of cos a; as well: 


cos (5 = — 


2 COS^ 023 


sin 2012 sin 2023 sin 0 i 3 


(1 — 2 COS^ 013 COS^ 023) 


cos^ 0^" 


12 


COS^ 023 


— sin^ 012 tan^ 023 


-|- (cos^ 013 sin^ 05(2 — cos^ 0 i 2 sin^ 013 ) -|- y/l — 2 cos^ 0 i 3 cos^ 023 cos lo sin 20 

cos 023 


(60) 


sin^ 0^3 


and eqs. (57) and (60) read: 

\Uri\ = I sin023 sin0i2 — sin0i3 cos $23 cos012®*'^| 


COS^ 023 


I COS 05^2 sin 0®3 — e cos 0®3 sin 0^2 sin 0231 , 

(62) 


COS^ 023 


COS b = ——-— 

sin 2012 sin 2023 sin 013 ^ 

013 COS^ 823 , . 2 nv , 1 ( 

- cos 20^2 ( 1-2“Z7-^23 + 1) ) + COS 2012 ( COS 

V cos ^23 ' 


^ 023 sin^ 013 - sin^ 023 ) 


(63) 


It follows from the results for cosd obtained for cosw / 0, eqs. (60) and (63), that in the 
case analysed in the present subsection one can obtain predictions for cos b only in theoretical 
models in which the value of the phase u is fixed by the model. 


The Cases of 


- (^23^ ^13^ ^12) Rotations 


We consider next a generalisation of the cases analysed in Sectionj^with the presence of a third 
rotation matrix in Uu arising from the neutrino sector, i.e., we employ the parametrisation of 
U given in eq. m Non-zero values of 033 are inspired by certain types of flavour symmetries 

= — 7 r /4 and = sin“^(l/3), for instance, 


(see, e.g., 


44-47 ). In the case of 0^2 = 


^23 


we have the so-called tri-permuting (TP) pattern, which was proposed and studied in 44 


In the statistical analysis of the predictions for cos <5, b and the Jcp factor we will perform 
in Section we will consider three representative values of 0^3 discussed in the literature: 
^13 = '?r/20, vr/lO and sin“^(l/3). 
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For the parametrisation of the matrix U given in eq. ( 18 ) with [ij) = ( 23 ), no constraints 
on the phase 6 can be obtained. Indeed, after we recast U in the form 


U = R23{023) Ql i?13(0r3) Ru{ 0’(2) Qo , 


( 64 ) 


where sin"^ 623 and Qi are given in eqs. ( 35 ) and ( 36 ), respectively, we find employing a similar 
procedure used in the previous sections: 


sin^ 013 = sin^ 0i3 , sin"^ 023 = sin^ 623 , sin^ 0i2 = sin"^ 0i2 , sin <5 = sin /? . 


2 


( 65 ) 


Thus, there is no correlation between the Dirac CPV phase 6 and the mixing angles in this 
set-up. 


4.1 The Scheme with 0^2 ~ (^ 23 ; ^ 13 ) ^ 12 ) Rotations 


In the parametrisation of the matrix U given in eq. ( 18 ) with (ij) = ( 12 ), the phase oj in the 
matrix T is unphysical (it “commutes” with i?i2(0f2) and can be absorbed by the field). 
Hence, the matrix T contains only one physical phase (/>,'!' = diag(l,e*‘^, 1 ), and 4> = 


Taking this into account and using eq. ( 18 ) with (ij) = ( 12 ) and 0^3 = —7r/4, we get the 
following expressions for sin^ 0i3, sin^ 023 and sin^ 0i2: 


1 


sim 013 = |C/e 3 r = 2 013 + COS^ 012 013 — -^12 sin 013 


sin^ 023 = 


sin^ 012 = 


\U^.3? 

l-\Ue3\^ 

\Ue2? 


= 1 - 


COS^ 013 


2 (1 - sin2 0i3) ’ 


- sin"^ 012 (cos 012 + sin 0i2 sin 0i3 


u \2 


l-|I7e3P l-sin^0i3[2 
+ cos^ 012 cos^ 013 sin^ 012 + -^12 sin 0i2 (cos 0i2 + sin 0i2 sin 013) 


( 66 ) 

( 67 ) 

( 68 ) 


\/2 

Solving eq. (66) for X12 and inserting the solution in eq. (68), 
of 013 , 032, 013 and 032: 

. asin2 0f2 + /3 

sm 012 = 

Here the parameters 


i — sm (713 

and f3 are given by: 


a = 


1 


2 cos 2012 + sin 20i2 


13 = sin 012 


cos^ 013 sin 012 


nu , ay (■ Sin2 0i3' 

012+cos 012 (^sm0i3-^^ 


we find sin^ 0i2 as a function 

( 70 ) 


(71) 

(72) 


Inverting the formula for sin^ 0i2 allows us to express sin^ 032 in terms of 0i2, 0i3, 032, 0i3: 

1 2d\^ sin^ 033 - sin 20^2 sin^ 0i3 7.^01 

r.r.^2av ■ 


sin^ 032 = 


2 cos^ 013 sin033(sin^ 0i2 — sin^ 0i2) + sin 2032 013 “ 

cos 2032 013 012 012 013 
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In the limit of vanishing 6 '(^ we have sin^ 0 f 2 = 2 sin^ 0 i 3 , which corresponds to the case of 
negligible considered in po]. 


Using eq. ( 68 ), one can express cos(j) in terms of the “standard” mixing angles 612 ■, ^13 


and the angles 0 f 2 , 9 12 and which are assumed to have known values: 


cos 0 = 


2 cos^ 013 (sin0^2) ^(sin0i2) ^ sin^ 0i2 — 2 cos^ cot^ 0®2 “ (cot 0^2 + sin ^ig)^ 


X(cos 0 


5 ( 3 ) ^tan 0 


12 


2-v/2(cot 05^2 + sin0 


1 -1 


13l 


(74) 


We note that from the requirements (0 < sin^ 0f2 < 1) A (—1 < cos </><!) one can obtain for 
a given each of the symmetry values of 05^2 considered and 023 = —7r/4, lowe r and upper 
bounds on the value of sin^ 0i2. These bounds will be discussed in subsection 5.2 Comparing 


the expressions for Jcp = lm(U*iU*^Ue 3 Ufj,i), obtained using eq. (18) with (ij) = (12) and 
023 = —vr/4, and in the standard parametrisation of U, one gets tne relation between sinc/i 
and sin 5: 


sin 5 = — 


sin 20?, 


12 


2 sin 2012 sin 20 i 3 sin 623 


cos^ 0?3 sin 20?2 + 2 cos 20?2 sin 0^3 


sm( 


(75) 


Similarly to the method employed in the previous Section, we use the equality of the ex¬ 
pressions for \Uri\ in the two parametrisations in order to derive the sum rule for cosJ of 
interest: 


\Uti\ = I sin023 sin 0 i 2 — sin 0 i 3 cos 623 cos 0 i 2 e*'^| = -j= \ sin 0?2 + cos 0?2 sin0?31 . (76) 

V2 


1 


From the above equation we find the following sum rule for cos <5: 

1 


cosd = 


(1 — 2 sin^ 013 + sin^ 05(3) sin^ 0i2 

(77) 


sin 2012 sin 0 i 3 | cos 0 ? 3 | (1 — 2 sin^ 0 i 3 -|- sin^ 0 ^ 3 ) 2 . 

+ cos^ 012 sin^ 013 cos^ 0?3 — cos^ 013 (sin 0?2 + cos 0?'2 sin 0 ? 3 )^ 


For 0?3 = 0 this sum rule reduces to the sum rule for cos<5 given in eq. (50) in plC)| . 

In the parametrisation of the PMNS matrix considered in this subsection, the rephasing 
invariant Jcp has the form: 


>7cp = — 


1 


8 V 2 


sin (j) cos 0^3 sin 20?2 [cos^ 0?3 sin 20^2 + 2 sin 0^3 cos 205(2] • (^ 8 ) 


In the case when 023 has a fixed value which differs from —vr/d, the expression for sin^ 623 , 

cos^ 023 cos^ 0?q 


eq. (67), changes as follows: 


sin 023 = 


\U^3\^ 

l-\Ue3\^ 


Equations (76) and (77) are also modified: 


= 1 - 


\Uri\ = I sin 023 sin 012 — sin 013 COS 023 COS 0126*^1 = I sin 0^2 sin 023 — cos 023 cos 05(2 sin 0 ^ 31 , 


As 


1 — sim 013 


(79) 


(80) 
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and 


X 


(cos^ 013 — cos^ 013 cos" 




- COS^ 013 (cos 012 sin 05(3 


0(^3 cos 023 ■ 


- sin 032 sin 023 )^ 


( 81 ) 


In the case of bi-trimaximal mixing 46 , i.e., for 0^2 = 023 = tan ^(V 3 — 1 ) and 0(^3 = 


48 


sin ^((3 — -v/ 3 )/ 6 ), the sum rule we have derived reduces to the sum rule obtained in 
However, this case is statistically disfavored by the current global neutrino oscillation data. 


4.2 The Scheme with 0^3 — (023)^13)^12) Rotations 

Here we switch to the parametrisation of the matrix U given in eq. ( | 18 | ) with [ij) 
Now the phase ijj in the matrix 'h is unphysical, and di = diag(l, 1 , Fixing 023 

and using also the standard parametrisation of U, we find: 

sin^ 013 = |t/g3|2 = 1 sin2 0^3 cos^ 0i3 + cos^ 013 sin^ 0i3 + X13 sin0^3 , 


sin^ 023 = 




COS^ 033 


l-|f7e3|2 2(1 -sin2 033) ’ 


sim 012 = 


\Ue2? 


X sin^ 013 (cos 012 - sin 032 sin 033)" 


l-|fle3P l-sin^0i3[2 
+ COS^ 033 COS^ 033 sin^ 032 + -^13 sin 032 (cos 032 — sin 032 sin 033) 


) ( 13 ). 
—7r/4 


( 82 ) 

( 83 ) 

( 84 ) 


Here 


Xi3 = —— sin 2033 *^0® ^13 ^ 

V2 

not dificult to find sin' 


Solving eq. ( 82 ) for X13 and inserting the solution in eq. ( 84 ), it is 
as a function of 0i3, 032, 033 and 033: 

. ps\v?ei^ + ri 

sm 012 = -r^-— , 

1 — sm^ 033 


where p and p are given by 

2 cos 2032 


1 

^=4 


COS^ 033 


sin 033 


V 


= sin 032 


*((2-sin20^2- 

cos^ 013 sin 032 - cos 032 ( sin 033 - 

V ®l'^ ^13 


sin^ 013 ' 


( 85 ) 

( 86 ) 

( 87 ) 

( 88 ) 


Using eq. (86) for sin^ 0i2 with p and p as given above, one can express sin^ 0f3 in terms of 
012, 013, 012) 6^13: 

. 2 QB _ 2 cos^ 013 sin 033(sin^ 012 — sin^ 05(2) — sin 2032 sin^ 013 + sin 2032 sin^ 033 /on' 

cos 2032 ®™ ^13 “ ^O® ^12 ®™ ^12 ^°®^ ^13 
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In the limit of vanishing 6 ^^ we find sin^ = 2 sin^ @13, as obtained in subsection 2.2 


Further, using eq. (84), we can write costa in terms of the standard parametrisation mixing 


angles and the known 0^2 and 61 ^: 


costa 


Icos^ 0i3(sin033) ^(sin032) ^ sin^ 0i2 — 2 cos^ 013 cot^ 0^3 — (cot0i2 “Sin0i3)^ 

1 -1 


X (cos 013) ^ tan 013 


2-v/2(cot 012 — sin 013) 


(90) 


Analogously to the case considered in the preceding subsection, from the requirements 
(0 < sin^ 0f3 < 1)A(—1 < costa < 1) one can obtain for a given 0i3, each of the symmetry 
values of 0i2 considered and 023 = — 7r/4 lower and upper bounds on the value of sin^0i2. 
These bounds will be discussed in subsection 15.31 

Comparing again the imaginary parts of U*iU*^Ue 3 Ufj,i, obtained using eq. (18) with 
{ij) = (13) and 9^^ = —vr/4, and in the standard parametrisation of U, one gets the following 
relation between sinta and sin 5 for arbitrarily fixed 0i2 and 013: 

sin 20f3 


sin S ^ __ 

■ 2012 sin 2013 cos 023 


2 sin 

Exploiting the equality of the expressions 


013 sin 2012 “ 2 cos 2012 sin 0i3 


cos 013 

for \Ufj_i\ written in 

1 


sinta. (91) 

the two parametrisations, 

U^i I = I cos 023 sin 012 + cos 0i2 sin 0i3 sin 623 \ = \ cos 0i2 sin 0i3 — sin 0i21 , (92) 

V2 


we get 


cos (5 = — 


the following sum rule for cos <5: 

1 


sin2012 sin0i3| cos0i3|(l — 2 sin^ 0i3 + sin^ ^13) ^ 

+ cos^ 012 sin^ 013 cos^ 0i3 — cos^ 013 (sin 0i2 — cos 0i2 sin 013)^ 


(1 — 2 sin^ 013 + sin^ 013) sin^ 0i2 

(93) 


For 013 = 0 this sum rule reduces to the sum rule for cos <5 given in eq. (27). 

In the parametrisation of the PMNS matrix considered in this subsection, the Jcp factor 
reads: 

Jcp =-^ sin ta cos 0i3 sin 20^3 [cos^ 0i3 sin 20i2 — 2 sin 0i3 cos 205(2] • (9“^) 

8v 2 

In the case of an arbitrary fixed value of 023, as it is not difficult to show, we have: 


and 


sin^ 023 = 


I = I cos 023 sin 012 + cos 0i2 sin 0i3 sin 0231 = | cos 0i2 sin 0i3 sin 623 + sin 012 cos 0231 


\U^.3? 


sin^ 023 cos^ 05(3 


l-|^e3p l-sin2 0i3 


(95) 


(96) 
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Using eqs. (95) and (96), we obtain in this case 


X 


"23J 


--J.JI- 

(cos^ 013 — cos^ 013 sin 
3^ 013 (cos 012 sin 05(3 sin 023 + sin 0^2 


- cos' 


0(^2 COS 023) 


( 97 ) 


The sum rules derived in Sections H-i and correspondin g to arbitrary hxed values of 
the angles contained in the matrix U^,, eqs. (13), (21), (32), (45), (63), (81) and (97), are 
summarised in Table In Table we give the corresponding formulae for sin^ 623 . 
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Parametrisation of U 

COS 5 

i? 12 ( 0 f 2 )^^ 23 ( 02 " 3 )^ 12 W 2 )Qo 

(cos 2013 — cos 2023) ^ 

COS 205^2 + (sin^ 012 — cos^ 012) 

2 sin^ 023 — (3 + cos 20^3) sin^ 0i3 


\/2 sin 2012 sin 0i3 cos 0^3 

cos 2013 ~ cos 2023 


(cos 2013 + cos 2023) ^ 

cos 205^2 + (sin^ 012 — cos^ 012) 

2 COS^ 023 — (3 — cos 2023 ) ^13 


\/2 sin 2012 sin 0i3 sin 0^3 

cos 2013 + cos 2023 

Rl 2 {ei 2 ) R 23 { 0 h) ^ R 23 { 0 h) Ri 2 { 9 ’( 2 ) Qo 

tan 023 

sin 2012 sin 0i3 

[cos 2032 + (sin^ 9 i 2 — cos^ 012) 

(1 - cot^ 023 sin^0l3)] 

Rl 3 { 9 h) R 23 { 9 h) ^ R 23 { 9 h) ^12(0^2) Qo 

cot 023 

sin 2012 sin 013 

[cos 2032 + (sin^ 012 — cos^ 032) 

(1 - tan^ 023 sin^ 013)] 

^12(012) ^ 13 ( 0 f 3 ) ^ R 23 { 9 ^ 23 ) Rl 2 { 9 'i^) Qo 

1 

2fi; cos w sin 20(^2 sin 023 cos 013 cos 023 , 9^,. 9„ 2„ 

3)^ 

sin 2012 sin 2023 sin 0i3 
/ cos^ 

— cos 20i^n ( 1- 

V C( 

2 nu (723 COb f 7 i 3 COb U 2 

COS ^23 

2 (sm^ 023 + 1 )) + cos 2012 (COS^ 023 sin^ 013 sin^ 023) 

3 S ^23 2 

Rl 2 { 9 h) ^ R 23 i 9 h) Ri 3 { 9 ’( 3 ) Ri 2 { 9 ’( 2 ) Qo 

1 

fl. „ r,r ^^2 QV QV \^\„2 

sin 2012 sin 013 cos 05(3 cos023 (cos^ 0i3 — cos^ 05^3 cos^ ^23)^ 
+ cos^ 012 sin^ 013 cos^ 0i3 cos^ 023 — cos^ 013 (cos 012 sin 

013 cos 023 — sin 032 sin 0^3)^ 

Olli Wl'l 

Rl 3 i 9 ts) ^ R 23 { 9 ’^ 23 ) ^13(0^3) ^ 12 (^r 2 ) Qo 

1 

r,r ^^2 QU ^\^2 QV \ ^-2 Q, ^ 

sin 2012 sin 013 cos 033 sin023 (cos^ 0i3 — cos^ 9 '(^ sin^ ^23)^ 
+ cos^ 012 sin^ 013 cos^ 05^3 sin^ 023 — cos^ 0i3(cos 05(2 sin 

033 sin 023 + sin 032 cos 023 )^ 

Dii.1 


Table 1: Summary of the sum rules for cosi5. The parameter k is defined in Section 3.2 after eq. (58 
parametrisation of U, i?i2(0i2).R23(023)^^23(023)-^i2(^i2)Qo! is the one quoted in eq. (13 and was derived in [30 


The sum rule corresponding to the 





































Parametrisation of U 

sin^ 023 

Rl2{0t2)-^R23m3)Rl2{0’(2)Q0 

sin^ 023 ~ sin^ 0 i 3 

1 — sin^ 013 

Rl3{Ot3)'^R23{0h)Rl2iO’(2)Qo 

sin^ 0^3 

1 — sin^ 013 

Rl2i9t2) R23{0h) ^ R23i9h) ^ 12 ( 0 ^ 2 ) QO 

sin^ 023 - sin^ 0 i 3 

1 — sin^ 013 

Rl3i9t3) R23{9h) ^ R23{9^23) Rl2{9'i^) Qo 

sin^ 023 

1 — sin^ 013 

Ri 2{9\2) Rim3) R23i9h) ^ 12 ( 0 ^ 2 ) Qo 

sin^ 023 — sin^ 0 i 3 + sin^ 0^3 cos^ 023 

1 — sin^ 013 

Rl2i9t2) ^ R23{9^23) RM9 i3) ^ 12 ( 0 ^ 2 ) Qo 

^ COS^ 023 COS^ 013 

1 — sin^ 013 

Rn{9t3) T R23{9-23) Risi^h) ^ 12 ( 0 ^ 2 ) Qo 

sin^ 023 COS^ 013 

1 — sin^ 013 


Table 2: Summary of the formulae for sin^ 6 * 23 . 


The formula for sin'^ 023 is given in eq. (35). 


5 Predictions 


In this Section we present results of a statistical analysis, performed using the procedure 
described in Appendix [a| (see also 34 ), which allows us to get the dependence of the 
function on the value of 6 and on the value of the Jcp factor. In what follows we always 


assume that ^23 = —vr/d. We find that in the case corresponding to eq. (16) with (ij) = (12), 
analysed in 30 , the results for as a function of 5 or Jcp are rather similar to those obtained 


in in the case of the parametrisation defined by eq. ( [ItI ) with (ij) — (kl) = (12) — (23). The 
main difference between these two cases is the predictions for sin^ 623 , which can deviate only 
by approximately 0.5 sin^ 813 from 0.5 in the first case and by a significantly larger amount in 
the second. As a consequence, the predictions in the first case are somewhat less favoured by 
the current data than in the second case, which is reflected in the higher value of ^it the 
Similar conclusions hold comparing the results in the case of 9f^ — { 623 , ^12) 
and in the corresponding case defined by eq. 0 


minimum. 


A mm 


rotations, described in Section 2.2 


with 

(ij) — [kl) = (13) — (23) and discussed in Section 3.1 Therefore, in what concerns these four 
schemes, in what follows we will present results of the statistical analysis of the predictions 
for 6 and the Jcp factor only for the scheme with (^fg rotations, considered 

in Section ISTI 
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We show in Tables and the predictions for cos <5 and 5 for all the schemes considered in 
the present study using the current best fit values of the neutrino mixing parameters sin^ 612 , 
sin^ 023 and sin^0i3, quoted in eqs. - (^, which enter into the sum rule expressions for 
cos (5, eqs. (13), (27), (45), (58), ([^), ( 1 ^) and eq. (50) in ref. 30 , unless other values of 


the indicated mixing parameters are explicitly specified. We present results only for the NO 
neutrino mass spectrum, since the results for the 10 spectrum differ insignificantly. Several 
comments are in order. 


Scheme 

TBM 

GRA 

GRB 

HG 

BM (LG) 


-0.114 

0.289 

-0.200 

0.476 

— 

^13 “ (^23’ ^ 12 ) 

0.114 

-0.289 

0.200 

-0.476 

— 

m2.oh)-{e^2M 

-0.091 

0.275 

-0.169 

0.445 

— 

(^f3)^i3)-(023)0r2) 

0.151 

-0.315 

0.251 

-0.531 

— 

(^!2)0f3)-(023)^r2) 

-0.122 

0.282 

-0.208 

0.469 

— 

Scheme 

[7r/20, —vr/4] 

[tt/IO, —vr/4] 

[a, -7r/4] 

[vr/20,6] 

[7r/20,7r/6] 

^12 “ (^23) ^13) ^ 12 ) 

-0.222 

0.760 

0.911 

-0.775 

-0.562 

Scheme 

[vr/20,c] 

[7r/20,7r/4] 

[vr/lO, 7r/4] 

[a,7r/4] 

[7r/20,d] 

^13 “ (^23) ^13) ^ 12 ) 

-0.866 

0.222 

-0.760 

-0.911 

-0.791 


Table 3: The predicted values of cos 5 using the current best fit values of the mixing angles, 
quoted in eqs. (§-(§ and corresponding to neutrino mass spectrum with NO, except for the 
case (012, ^is) — (023) ^ 12 ) whh w = 0 and k = 1 , in which sin^ 023 = 0.48802 is used. We have 
defined a = sin“^(l/3), b = sin“^(l/-\/2 + r), c = sin“^(l/-\/3) and d = sin“^(V3 — r/ 2 ). For 
the last two schemes we give in square brackets the values of [0i3,0i2]- TBM, GRA, GRB, 
HG and BM (LG) refer, in particular, to the different fixed values of 0i2 = c, 6, d, tt/O and 7r/4, 
respectively. See text for further details. 


Scheme 

TBM 

GRA 

GRB 

HG 

BM (LG) 

^12 “ (^23’ ^12) 

97 V 263 

73 V 287 

102 V 258 

62 V 298 

— 

Oh-{0h:(^12) 

83 V 277 

107 V 253 

78 V 282 

118 V 242 

— 

i0t2,0h)-i0h,012) 

95 V 265 

74 V 286 

100 V 260 

64 V 296 

— 

i0t3,0h)-i0h,012) 

81 V 279 

108 V 252 

75 V 285 

122 V 238 

— 

idt2,0t3)-i0h,012) 

97 V 263 

74 V 286 

102 V 258 

62 V 298 

— 

Scheme 

[7r/20, —vr/4] 

[vr/10, —vr/4] 

[a, -vr/4] 

[vr/20, 6] 

[vr/20, vr/6] 

^12 “ 1^23’ ^135 ^ 12 ) 

103 V 257 

41 V 319 

24 V 336 

141 V 219 

124 V 236 

Scheme 

[vr/20,c] 

[vr/20, vr/4] 

[vr/10, vr/4] 

[a, vr/4] 

[vr/20, d] 

oe /nu r\v f\v \ 

^13 “ 1^23’ ^135 ^ 12 ) 

150 V 210 

77 V 283 

139 V 221 

156 V 204 

142 V 218 


Table 4: The same as in Tablebut for 6 given in degrees (see text for further details). 


We do not present predictions for the BM (LC) symmetry form of 11^ in Tables and 
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because for the current best fit values of sin^ 0i2, sin 
give unphysical values of cos 5 (see, however, refs. 


^ 023 , sin^ 013 the corresponding sum rules 
Using the best fit value of sin^ 0i3, 


30 


34 




^12), 


we get physical values of cos 6 in the BM case for the following minimal values of sin^ 0i2: 
cos 5 = —0.993 (5 = tt) for sin^ 0 i 2 = 0.348 in the scheme 6I2 — (^23)^12)1 

cos 5 = +0.993 (5 + 0) for sin^ 0i2 = 0.348 in the scheme 0f3 — {O^s, ^ 12 )) 

cos 5 = —0.994 (5 + tt) for sin^ 0i2 = 0.349 in the scheme { 0 I 21 ^ 13 ) “ (^ 23 ;' 

cos 5 = —0.996 (5 + tt) for sin^ 0i2 = 0.332 in the scheme (0^2) ^ 23 ) “ (^ 23 ;' 

cos 5 = +0.997 {5 = 0) for sin^ 0i2 = 0.368 in the scheme ( 0 ^ 3 , 013 ) — (02 

where in the case of the scheme (0f2, ^13) (^23 >^12) we fixed sin^ 623 = 0.48802 (we will 
comment later on this choice), while sin^ 623 was set to its best fit value for the last two 
set-ups. 

Results for the scheme (012,^23) “ (^23) in the cases of the TBM and BM symmetry 
forms of the matrix Ui, were presented first in 
the GRA, GRB and HG symmetry forms of U^, 


+3> ^12) 


31 , while results for the same scheme and 
as well as for the scheme 0^2 — (^23) ^12) 

The predictions for cos <5 and 6 


30 


all symmetry forms considered, were obtained hrst in 
were derived in 30 and for the best fit values of the relevant neutrino mixing parameters 
found in an earlier global analysis performed in [3^ and differ somewhat (albeit not much) 


from those quoted in Tables and The values under discussion given in these tables are 
from and correspond to the best fit values quoted in eqs. (§-(§. 

The predictions for cos <5 of the 0®2 — (023; ^12) ^13 “ (^23> ^12) schemes for each of the 

symmetry forms of Uiy considered differ only by sign. The 0f2 — (^235^12) scheme and the 
(^12’ ^13) “ (^23’ ^12) scheme with a; = 0 provide very similar predictions for cos 5. 

In the schemes with three rotations in we consider, cosJ has values which differ sig- 
nihcantly (being larger in absolute value) from the values predicted by the schemes with two 
rotations in discussed by us, the only exceptions being i) the 0^2(13) “ i^ 23 i ^13) ^12) scheme 
with [013,012] = 1^/20,^/4], for which |cos5| = 0.222, and ii) 0i2 — (^23)^13)^12) scheme 
with [013,012] = [vr/20,7r/6] in which cos 5 = —0.562. 

The predictions for cos (5 of the schemes denoted as (012,^23) - (^23) ^12 ) and (0f3,013) - 
(023, 012) differ for each of the symmetry forms of Ui, considered both by sign and magnitude. 
If the best fit value of 023 were 7r/4, these predictions would differ only by sign. 

In the case of the (012,^13) -(^2^3) ^1^2) scheme with a; = 0, the predictions for cosd are 
very sensitive to the value of sin^ 023. Using the best fit values of sin^ 0i2 and sin^ 0i3 for 
the NO neutrino mass spectrum, quoted in eqs. (© and ([^, we find from the constraints 
( —1 < cos'i/' < 1) and (0 < sin^ 0f3 < 1) A (0 < sin^ 0f2 < 1), where sin^ 0^3, sin^ 0f2 and cosi/) 


are given in eqs. (53) - (55), that sim 023 should lie in the following intervals: 


(0.488,0.496) U (0.847,0.909) for TBM; 

(0.488,0.519) U (0.948,0.971) for BM; 

(0.488,0.497) U (0.807,0.880) for GRA; 

(0.488,0.498) U (0.856,0.914) for GRB; 

(0.488,0.500) U (0.787,0.866) for HG. 

Obviously, the quoted intervals with sin^ 023 > 0.78 are ruled out by the current data. We 
observe that a small increase of sin^ 023 from the value 0.48802 produces a relatively large 

®For sin^ O23 < 0.48802, cos 5 has an unphysical (complex) value. 
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variation of cos 6 . The strong dependence of cos 6 on sin^ 623 takes place for values of io 
satisfying roughly cos a; ^ 0.01. In contrast, for cos to = 0, cos exhibits a relatively weak 
dependence on sin^ 623 ■ For the reasons related to the dependence of cos <5 on a; we are not 
going to present results of the statistical analysis in this case. This can be done in specific 
models of neutrino mixing, in which the value of the phase oj is fixed by the model. 


5.1 The Scheme with (^13,^23) ~ (^23)^12) Rotations 
In the left panel of Fig. we show the likelihood function, defined as 

Tf ( X^(cos(i)' 

L(cosoj oc exp- 


(98) 


versus cos <5 for the NO neutrino mass spectrum for the scheme with (013,^23) - (^23, ^12) 
This function represents the most probable values of cos 5 for each of the sym- 


10 


rotations 

metry forms considered. In the analysis performed by us we use as input the current global 
neutrino oscillation data on sin^0i2, sin^023) sin^ @13 and <5 15 . The maxima of L{cos 6 ), 
= Xmin)’ different symmetry forms of Uu considered, correspond to the values 

of cos 6 given in Table The results shown are obtained by marginalising over sin^ 0i3 and 


34 


sin^ 023 for a fixed value of S (for details of the statistical analysis see Appendix and 
The na confidence level (C.L.) region corresponds to the interval of values of cos 5 for which 
L(cosS) > L(x^ = Xmin) ‘ Here Xmin value of x^ in the minimum. 

As can be observed from the left panel of Fig. for the TBM and GRB forms there is a 
substantial overlap of the corresponding likelihood functions. The same observation holds also 
for the GRA and HG forms. However, the likelihood functions of these two sets of symmetry 
forms overlap only at Sa and in a small interval of values of cos <5. Thus, the TBM/GRB, 
GRA/HG and BM (LG) symmetry forms might be distinguished with a not very demanding 
(in terms of precision) measurement of cos <5. At the maximum, the non-normalised likelihood 
function equals exp(—x^in/^); this value allows one to judge quantitatively about the 
compatibility of a given symmetry form with the global neutrino oscillation data, as we have 
pointed out. 

In the right panel of Fig. we present L versus cos <5 within the Gaussian approximation 
for details), using the current best fit values of sin^0i2, sin^023, sin^ ^13 for NO 
spectrum, given in eqs. (©-(©, and the prospective la uncertainties in the measurement 
of these mixing parameters. More specifically, we use as Icr uncertainties i) 0.7% for sin^ ^12, 
which is the prospective sensitivity of the JUNO experiment 49 , ii) 5% for sin^ @23 


(see 34 


obtained from the prospective uncertainty of 2% on sin^ 2023 expected to be reached in 
the NOvA and T2K experiments, and iii) 3% for sin^ ^13, deduced from the error of 3% 
on sin^ 2013 planned to be reached in the Daya Bay experiment [3,51 . The BM (LG) case 
is quite sensitive to the values of sin^ 0i2 and sin^ 023 and for the current best fit values is 
disfavoured at more than 2a. That the BM (LG) case is disfavoured by the current data 
can be understood, in particular, from the following observation. Using the best fit values 
of sin^ 013 and sin^ 0i2 as well as the constraint — 1 < cos a < 1 , where cos a is defined in 
eq. (42), one finds that sin^ 023 should satisfy sin^ 023 > 0.63, which practically coincides with 
the currently allowed maximal value of sin^ 023 at 3a (see eq. (0). 

^®The corresponding results for the IO neutrino mass spectrum differ little from those shown in the left 
panel of Fig. 


^This sensitivity is planned to be achieved in future neutrino facilities 
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Figure 1: The likelihood function versus cos 5 for the NO neutrino mass spectrum after 
marginalising over sin^ @13 and sin^ 623 for the TBM, BM (LC), GRA, GRB and HG symmetry 
forms of the matrix in the (^ 13 , 6 * 23 ) ~ (^ 23 ; ^ 12 ) set-up. The results shown are obtained 
using eq. (451 and i) the latest results on the mixing parameters sin^ 612 , sin^ ^ 13 , sin^ 623 and 
6 found in the global analysis of the neutrino oscillation data (left panel), and ii) the 
prospective Icr uncertainties on sin^ 6*12, sin^ 613, sin^ 623 and the Gaussian approximation for 


the likelihood function (right panel) (see text for further details). 


It is interesting to compare the results described above and obtained in the scheme denoted 
by (013,^23) “ (^23)^12) those obtained in in the (^12, ^23) “ (^235^12) set-up. We 
recall that for each of the symmetry forms we have considered — TBM, BM, GRA, GRB 
and HG — 9\2 has a specific fixed value and ^23 = —7r/4. The first thing to note is that for 
a given symmetry form, cos <5 is predicted to have opposite signs in the two schemes. In the 
scheme ( 013 , ^Ig) - analysed in the present article, one has cos <5 > 0 in the TBM, 

GRB and BM (LC) cases, while cos <5 < 0 in the cases of the GRA and HG symmetry forms. 
As in the (^12,^23) “ (^23)^12) set-up, there are significant overlaps between the TBM/GRB 
and GRA/HG forms of respectively. The BM (LC) case is disfavoured at more than 2(j 
confidence level. It is also important to notice that due to the fact that the best fit value of 
sin^ 023 < 0.5, the predictions for cos <5 for each symmetry form, obtained in the two set-ups 


differ not only by sign but also in absolute value, as was already pointed out in Section 3.1 


Thus, a precise measurement of cos 5 would allow one to distinguish not only between the 
symmetry forms of but also could provide an indication about the structure of the matrix 
Ue.. 


We note that the predictions for sin^ 023 are rather similar in the cases of the two schemes 
discussed, (0f3, 023) — (023, 012) and (0^, 023) ~ (^235 ^12)- completeness Nfj = 

as a function of sin^ 023 in Appendix 

For the rephasing invariant Jcp, using the current global neutrino oscillation data, we 
find for the symmetry forms considered the following best fit values and the 3ct ranges for the 
NO neutrino mass spectrum: 


-0.039 < Jcp < -0.026, 0.030 < Jcp < 0.036 for TBM; (99) 

-0.026 < Jcp < 0.023 for BM (LC) ; (100) 

-0.037 < Jcp < -0.024, 0.029 < Jcp < 0.035 for GRA; (101) 


Jcp = — 0.033 
Jcp = — 0.004 
Jcp = — 0.032 
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Jcp = - 0.033 , -0.039 < Jcp < -0.023, 0.028 < Jcp < 0.036 for GRB; (102) 

Jcp = - 0.028 , -0.035 < Jcp < -0.014, 0.021 < Jcp < 0.032 for HG. (103) 

Thus, relatively large GP-violating effects in neutrino oscillations are predicted for all sym¬ 
metry forms considered, the only exception being the case of the BM symmetry form. 


5.2 The Scheme with d\2 — (^23)^135^12) Rotations 

For the scheme with 612 — (^23)^13)^12) rotations we hnd that only for particular values of 
9 i 2 ^-rid among those considered by us, the allowed intervals of values of sin^ 612 satisfy 
the requirement that they contain in addition to the best fit value of siv? 612 also the 1.5a 
experimentally allowed range of sin^0i2. Indeed, combining the conditions 0 < sit? 6\2 < 1 
and |cos(/>| < 1, where sin^ 0f2 rmd cos4> are given in eqs. (73) and (74), respectively, and 


allowing sin^ ^13 to vary in the Ha range for the NO spectrum, we get restrictions on the value 
of sin^ 012, presented in Table We see from the Table that only five out of 18 combinations 
of the angles 0^2 ^rid 0^3 considered by us satisfy the requirement formulated above. In Table[^ 
these cases are marked with the subscripts I, II, III, IV, V, while the ones marked with an 
asterisk contain values of sin^ 0i2 allowed at 2a 15 . 


0’(2 

0^3 = vr /20 

0^3 = vr/lO 

^13 = sin ^1/3) 

sin“^(l/-v/3) 

(0.319,0.654)* 

(0.471,0.773) 

(0.495,0.789) 

7 r /4 

(0.484,0.803) 

(0.639,0.897) 

(0.662,0.909) 

— 7 r /4 

(0.197,0.516)iii 

(0.103,0.361)1 

(0.091,0.338)iv 

sin“^(l/-v/2 + r) 

(0.262,0.594)11 

(0.409,0.719) 

(0.434,0.737) 

sin“^(\/3 — r/2) 

(0.331,0.666)* 

(0.484,0.784) 

(0.508,0.800) 

7 r /6 

(0.236,0.564)v 

(0.380,0.692) 

(0.404,0.710) 


Table 5: Ranges of sin^ 0i2 obtained from the requirements (0 < sin^ 0®2 < 1) ^ (~1 < cos 4> < 
1) allowing sin^ 0i3 to vary in the 3cr allowed range for the NO neutrino mass spectrum, quoted 
in eq. The cases for which the best fit value of sin^ 0i2 = 0.308 is within the corresponding 
allowed ranges are marked with the subscripts I, II, III, IV, V. The cases marked with an 
asterisk contain values of sin^ 0i2 allowed at 2a 15 . 


Equation (67) implies that sin^ 023 is fixed by the value of 0i3, and for the best fit value 
of sin^ 013 and the values of 0^3 = 0, 7r/20 , tt/IO, sin“^(l/3), considered by us, we get, 
respectively: sin^ 023 = 0.488, 0.501, 0.537, 0.545. Therefore a measurement of sin^ 023 with 
a sufficiently high precision would rule out at least some of the cases with fixed values of 9'^^ 
considered in the literature. 

We will perform a statistical analysis of the predictions for cos 6 in the five cases — I, II, 
III, IV, V — listed above. The analysis is similar to the one discussed in Section 5.1 The 


only difference is that when we consider the prospective sensitivities on the PMNS mixing 
angles we will assume sin^ 023 to have the following potential best fit values: sin^ 023 = 0.488, 
0.501, 0.537, 0.545. Note that for the best fit value of sin^0i3, sin^ 023 = 0.488 does not 
correspond to any of the values of 0^3 in the five cases — I, II, III, IV, V — of interest. Thus, 
sin^ 023 = 0.488 is not the most probable value in any of the hve cases considered: depending 
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on the case, the most probable value is one of the other three values of sin^ 623 listed above. 
We include results for sin^ 623 = 0.488 to illustrate how the likelihood function changes when 
the best fit value of sin^ O23, determined in a global analysis, differs from the value of sin^ 623 
predicted in a given case. 

In Fig. [^we show the likelihood function versus cosd for all the cases marked with the 
subscripts in Table The maxima of the likelihood function in the hve cases considered take 
place at the corresponding values of cos 5 cited in Table As Fig. clearly indicates, the 
cases differ not only in the predictions for sin^ 623 , which in the considered set-up is a function 
of sin^ 05^3 and sin^ 613 , but also in the predictions for cosJ. Given the values of 612 and 613 , 
the positions of the peaks are determined by the values of 612 and 0^3. 

The Cases I and IV are disfavoured by the current data because the corresponding values 
of sin^ 023 = 0.537 and 0.545 are disfavoured. The Cases II, III and V are less favoured for the 
NO neutrino mass spectrum than for the 10 spectrum since sin^ 623 = 0.501 is less favoured 
for the first than for the second spectrum. 

In Fig. we show the predictions for cos <5 using the prospective precision in the measure¬ 
ment of sin^ 012, sin^ 0i3, sin^ 023, the best fit values for sin^ 0i2 and sin^ 0i3 as in eqs. Q and 
([^ and the potential best fit values of sin^ 023 = 0.488, 0.501, 0.537, 0.545. The values of 
sin^ 023 correspond in the scheme discussed to the best fit value of sin^ 0i3 in the cases which 
are compatible with the current 1.5cr range of allowed values of sin^ 0i2. The position of the 
peaks, obviously, does not depend explicitly on the assumed experimentally determined best 
fit value of sin^ 023. For the best fit value of sin^ 0i3 used, the corresponding sum rule for 
cos S depends on the given fixed value of 0 ^^, and via it, on the predicted value of sin^ 023 (see 
eqs. (67) and ([TT])). Therefore, the compatibility of a given case with the considered hypo¬ 
thetical data on sin^ 023 clearly depends on the assumed best fit value of sin^ 023 determined 
from the data. 



Figure 2: The likelihood function versus cos <5 for the NO (10) neutrino mass spectrum in 
the left (right) panel after marginalising over sin^ 0 i 3 for the scheme 0^2 ~ ( 023 j ^ 12 ) with 

[ 05 ( 3 , 012 ] fixed as [tt/IO, —7r/4] (Case I), [7r/20, 6] (Case II), [7r/20, —7r/4] (Case III), [a, —7r/4] 
(Case IV), [7r/20,7r/6] (Case V), where a = sin“^(l/3) and b = sin“^(I/-y2 -|- r), r being the 
golden ratio. The figure is obtained using the sum rule in eq. and the latest results on 
sin^ 012, sin^ 0i3, sin^ 623 and 6 from the global analysis of the neutrino oscillation data 15 . 


As the results shown in Fig. [^indicate, distinguishing between the Cases I/IV and the 
other three cases would not require exceedingly high precision measurement of cos 6 . Distin- 
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Figure 3: The likelihood function versus cos 6 for the NO neutrino mass spectrum in the same 
cases as in Fig. but using the Gaussian approximation with the prospective uncertainties in 

““ o o o o o 

the measurement of sin 612, sin ^ 13 , sin 023 , the best fit values for sin 0 i 2 and sin @13 as in 
eqs. ^ and ^ and the potential best fit values of sin^ 623 = 0.488, 0.501, 0.537, 0.545. Upper 
left (right) panel: sin^ 023 = 0.488 (0.501); lower left (right) panel: sin^ 023 = 0.537 ( 0.545). 


guishing between the Cases II, III and V would be more challenging in terms of the requisite 
precision on cos 5. In both cases the precision required will depend, in particular, on the ex¬ 
perimentally determined best fit value of cos 5. As Fig. [^also indicates, one of the discussed 
two groups of Cases might be strongly disfavoured by the best fit value of sin^ 023 determined 
in the future high precision experiments. 

We have performed also a statistical analysis of the predictions for the rephasing invariant 
Jcp, minimising for fixed values of Jqp- We give No- = -\/i? as a function of Jcp in Fig.|^ 
The dashed lines represent the results of the global fit [im, while the solid ones represent 


the results we obtain for each of the considered cases, minimising the value of in 6^2 for 
a fixed value of Jcp using eq. (78). The blue lines correspond to the NO neutrino mass 
spectrum, while the red ones are for the 10 spectrum. The value of x^ fo the minimum, 
which corresponds to the best fit value of Jcp predicted in the model, allows one to conclude 
about compatibility of this model with the global neutrino oscillation data. As it can be 
observed from Fig. the zero value of Jcp in the Cases III and V is excluded at more than 
3cr with respect to the confidence level of the corresponding minimum. Although in the other 
three cases the best fit values of Jcp are relatively large, as their numerical values quoted 
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below show, Jcp = 0 is only weakly disfavoured statistically. The best fit values and the 3a 
ranges of the rephasing invariant JcP) obtained for the NO neutrino mass spectrum using the 
current global neutrino oscillation data, in the five cases considered by us are given by: 


Jcp = “ 0.023 , 
Jcp = ~ 0.022 , 
Jcp = — 0.033 , 
Jcp = - 0.016 , 
Jcp = ~ 0.028 , 


-0.032 < Jcp < 0.029 for Case I; 

-0.035 < Jcp < 0.031 for Case II; 

-0.039 < Jcp < -0.025, 0.030 < Jcp < 0.036 for Case III; 
-0.028 < Jcp < 0.026 for Case IV; 

-0.037 < Jcp < -0.010, 0.018 < Jcp < 0.034 for Case V. 


(104) 

(105) 

(106) 

(107) 

(108) 




Figure 4: N„- = as a function of Jqp in the scheme 9^2 — (9^3, 9^^, 9 ^ 2 ) with 9 ^ 2 ] 
fixed as [tt/IO, —7r/4] (Case 1), [7r/20,6] (Case II), [7r/20, — tt/I] (Case III), [a, —7r/4] (Case 
IV), [7r/20,7r/6] (Case V), where a = sin“^(l/3) and b = sin“^(l/V2 + r), r being the golden 
ratio. The dashed lines represent the results of the global fit [^, while the solid ones represent 
the results we obtain in our set-up. The blue (red) lines are for the NO (10) neutrino mass 
spectrum. 
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5.3 The Scheme with — (023; ^ 13 , ^u) Rotations 


As in the set-up discussed in the subsection 5.2, we find for the scheme with — { 623 , 013 , 012 ) 
rotations that only particular values of 05^2 allow one to obtain the current best fit 

value of sin^ 0i2. Combining the requirements 0 < sin^ < 1 and | cosa;| < 1, where sin^ 
and cosw are given in eqs. (89) and (90), respectively, and allowing sin^ 0i3 to vary in its 3 ct 
allowed range corresponding to the NO spectrum, we get restrictions on the value of sin^ 0i2, 
presented in Table It follows from the results in Table that only for five out of 18 
combinations of the angles 0^2 and 61 ^, the best fit value of sin^ 0i2 = 0.308 and the 1.5a 
experimentally allowed interval of values of sin^ 0i2 are inside the allowed ranges. In Table 
these cases are marked with the subscripts I, II, III, IV, V, while in the case marked with an 
asterisk, the allowed range contains values of sin^ 0 i 2 allowed at 2cr 


15 


0^2 

II 

to 

0 

0^3 = ^/lO 

013 = sin HV3) 

sin“^(l/-v/3) 

(0.081, 0.348)iii 

(0.024,0.209) 

(0.019,0.189) 

7r/4 

(0.197,0.516)1 

(0.103,0.361)iv 

(0.091, 0.338)ii 

—7r/4 

(0.484,0.803) 

(0.639,0.897) 

(0.662,0.909) 

sin“-'-(l/-v /2 -b r) 

(0.051,0.291)* 

(0.009,0.161) 

(0.006,0.143) 

sin“^(-v/3 — r/ 2 ) 

(0.089, 0.361)v 

(0.028,0.220) 

( 0 . 022 , 0 . 200 ) 

7r/6 

(0.038,0.264) 

(0.004,0.140) 

(0.002,0.123) 


Table 6: Ranges of sin^ 612 obtained from the requirements (0 < sin^ < 1) A (—1 < cos w < 
1) allowing sin^ 613 to vary in the 3cr allowed range for the NO neutrino mass spectrum, quoted 
in eq. The cases for which the best fit value of sin^ 612 = 0.308 is within the corresponding 
allowed ranges are marked with the subscripts I, II, III, IV, V. The case marked with an 
asterisk contains values of sin^ 612 allowed at 2cr 15 . 


The values of sin^ 623 in this model depend on the reactor angle 0 i 3 and 0^3 through 
eq. (83). Using the best fit value of sin^ 0 i 3 for the NO spectrum and eq. (83), we find 
sin^023 = 0.512, 0.499, 0.463, 0.455 for 6 '(^ = 0, 7 r/ 20 , vr/lO, sin“^(l/3), respectively. Thus, 
in the scheme under discussion sin^ 623 decreases with the increase of 0 ^ 3 , which is in contrast 
to the behaviour of sin^ 623 in the set-up discussed in the preceding subsection. As we have 
already remarked, a measurement of sin^ 623 with a sufficiently high precision, or at least 
the determination of the octant of 023 , would allow one to exclude some of the values of 9i^ 
considered in the literature. 

The statistical analyses for 5 and Jcp performed in the present subsection are similar to 
those performed in the previous subsections. In particular, we show in Fig. the dependence 
of the likelihood function on cos (5 using the current knowledge on the PMNS mixing angles 
and the Dirac CPV phase from the latest global fit results. Due to the very narrow prediction 
for sin^ 623 in this set-up, the prospective sensitivity likelihood curve depends strongly on the 
assumed best fit value of sin^ 623 . For this reason we present in Fig. the predictions for cos 6 
using the prospective sensitivities on the mixing angles, the best fit values for sin^ 0 i 2 and 
sin^ 013 as in eqs. ([^ and ([^ and the potential best fit values of sin^ 023 = 0.512, 0.499, 0.463, 
0.455. We use the value of sin^ 023 = 0.512, corresponding to 0 i 3 = 0, for the same reason we 
used the value of sin^ 023 = 0.488 in the analysis in the preceding subsection, where we gave 
also a detailed explanation. 
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Figure 5: The likelihood function versus cos <5 for the NO (10) neutrino mass spectrum 
in the left (right) panel after marginalising over sin^ 9i^ for the scheme — (^23 > ^131^ 12 ) 
with [013,^12] fixed as [7r/20,7r/4] (Case I), [a, 7r/4] (Case II), [7r/20,c] (Case III), [tt/IO, 7 r/4] 


(Case IV), [7r/20,(i] (Case V). We have defined a = sin“^(l/3), c = sin“^(l/-\/3) and d = 
sin“^(-\/3 — r/2), r being the golden ratio. The figure is obtained using the sum rule in eq. (931 
and the latest results on sin^ 612 , sin^ di 3 , sin^ 023 and 5 from the global analysis of the neutrino 
oscillation data 


15 


As Fig. clearly shows, the position of the peaks does not depend on the assumed best 
fit value of sin^023- However, the height of the peaks reflects to what degree the model is 
disfavoured due to the difference between the assumed best fit value of sin^ 023 and the value 
predicted in the corresponding set-up. 

The results shown in Fig. [^clearly indicate that i) the measurement of cos(5 can allow 
one to distinguish between the Case I and the other four cases; ii) distinguishing between 
the Cases II/III and the Cases IV/V might be possible, but is very challenging in terms of 
the precision on cosd required to achieve that; and iii) distinguishing between the Cases II 
and III (the Cases IV and V) seems practically impossible. Some of, or even all, these cases 
would be strongly disfavoured if the best fit value of sin^ 023 determined with the assumed 
high precision in the future experiments were relatively large, say, sin^ 023 0.54. 

The results on the predictions for the rephasing invariant Jcp are presented in Fig. 
where we show the dependence of on Jcp- It follows from the results presented in 

Fig.0 in particular, that Jcp = 0 is excluded at more than Su with respect to the confidence 
level of the corresponding minimum only in the Case I. For the rephasing invariant Jcp, 
using the current global neutrino oscillation data, we find for the different cases considered 
the following best fit values and 3 cj ranges for the NO neutrino mass spectrum: 


Jcp = 

-0.033, 

-0.039 

< 

-^CP 

< 

-0.025, C 

).029 < Jcp < 0.037 for Case I; 

(109) 

Jcp = 

-0.016, 

-0.028 

< 

Jcp 

< 

0.025 for 

Case II; 

(110) 

Jcp = 

-0.018, 

-0.029 

< 

Jcp 

< 

0.026 for 

Case III; 

(111) 

Jcp = 

-0.023, 

-0.031 

< 

Jcp 

< 

0.029 for 

Case IV; 

(112) 

Jcp = 

-0.022, 

-0.030 

< 

Jcp 

< 

0.028 for 

Case V. 

(113) 
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Figure 6: The likelihood function versus cosi5 for the NO neutrino mass spectrum in the 
cases described in Fig.[^ but within the Gaussian approximation. The upper left (right) panel 
corresponds to the potential best fit value of sin^ 023 = 0.5 1 2 (0.499), while the lower left 
(right) panel is obtained for the potential best fit value of sin^ 023 = 0.463 (0.455); the best 
fit values of sin^ 0i2 and sin^ 0i3 correspond to those quoted in eqs. ^ and The figure is 
obtained using the prospective uncertainties in the values of sin^ 0 i 2 , sin^ 0 i 3 and sin^ 623 ■ 


6 Summary and Conclusions 

In the present article we have derived predictions for the Dirac phase 5 present in the 3x3 
unitary neutrino mixing matrix U = ulUy = {Ue)'^ '^Uy Qq, where Ue {Ue) and Uy {Uy) are 
3x3 unitary (CKM-like) matrices which arise from the diagonalisation, respectively, of the 
charged lepton and the neutrino mass matrices, and 'h and Qq are diagonal phase matrices 
each containing in the general case two physical CPV phases. The phases in the matrix 
Qo contribute to the Majorana phases in the PMNS matrix. After performing a systematic 
search, we have considered forms of Ue and Uy allowing us to express <5 as a function of the 
PMNS mixing angles, 0 i 2 , ^13 and 023) present in U, and the angles contained in Uy. We have 
derived such sum rules for cos <5 in the cases of forms for which the sum rules of interest do 
not exist in the literature. More specifically, we have derived new sum rules for cos <5 in the 
following cases: 

i) U = Ru{eh)'^R23{e^3)Ri2{0’(2)Qo (^fa - (023>^12) scheme), 

ii) U = i?i 3 ( 0 f 3 )^i? 23 ( 02 " 3 )^i 2 ( 0 r 2 )Qo Scheme), 
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Figure 7 : The same as in Fig. but for the scheme — (023, ^* 13 , ^ 12 ) with O 12 } given 
by [ 7 r/ 20 , 7 r/ 4 ] (Case I), [a, 7 r/ 4 ] (Case II), [ 7 r/ 20 ,c] (Case III), [ 7 r/ 10 , 7 r/ 4 ] (Case IV), [ 7 r/ 20 ,d] 
(Case V), where a = sin“^(I/3), c = sin“^(l/-\/3) and d = sin“^(-v/3 — r/2), r being the golden 
ratio. The dashed lines represent the results of the global fit [^, while the solid ones represent 
the results we obtain in our set-up. The blue (red) lines are for the NO (10) neutrino mass 
spectrum. 


ill) U = Ris{et,)R2M,)^R23ie-,^)Ri2{0’(2)Qo - (^ 23 >^r 2 ) Scheme), 

iv) U = Ri2{0t2)Ri3{0h)'^R23{9h)Ri2i0’(2)Qo " (^23>^r2) Scheme), 

v) U = Ri2{9t2)^R23{e^23)RiM3)Ri2i0i2)Qo " (^23> ^^3-^^2) Scheme), and 

Vi) U = Ri3i9h)^R23{Oh)Ri3{Oh)Ri2i9’(2)Qo {9h - {923,91^,612) scheme), 

where Rij are real orthogonal matrices describing rotations in the i-j plane, and and 
stand for the rotation angles contained in Uf. and Uy, respectively. In the sum rules 
cos (5 is expressed, in general, in terms of the three angles of the PMNS matrix, 612 , ^13 and 
023 , measured, e.g., in the neutrino oscillation experiments, and the angles in Uy, which are 
assumed to have fixed known values. In the case of the scheme iv), cos <5 depends in addition 
on an a priori unknown phase w, whose value can only be fixed in a self-consistent model of 
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neutrino mass generation. A summary of the sum rules derived in the present article is given 
in Table [TJ 

To obtain predictions for cos <5, 5 and the Jcp factor, which controls the magnitude of the 
CP-violating effects in neutrino oscillations, we have considered several forms of determined 
by, or associated with, symmetries, for which the angles in have specific values. More 
concretely, in the cases i) - iv), we have performed analyses for the TBM, BM (LC), GRA, 
GRB, and HG forms of U,^. For all these forms we have @23 = — vr/4 and 0^3 = 0. The 
forms differ by the value of the angle 612 , which for the different forms of interest was given 
in the Introduction. In the schemes v) and vi) with non-zero fixed values of 0 ^ 3 , which are 
also inspired by certain types of flavour symmetries, we have considered three representative 
values of discussed in the literature, 05^3 = 7r/20, vr/lO and a = sin“^(l/3), in combination 
with specific values of 05^2 — altogether five sets of different pairs of values of [ 61 ^, 9 '( 2 ] in each 
of the two schemes. They are given in Table 

We first obtained predictions for cos <5 and <5 using the current best fit values of sin^ 0 i 2 , 
sin^ 013 and sin^ 023 , given in eqs. @-§. They are summarised in Tables and The 
quoted values of cos<5 and 6 for the scheme iv) are for a; = 0. For completeness, in Tables]^ 
and 1 ^ we have presented results also for 

vii) the ( 0^2 — (^ 23 >^ 12 ) scheme (in which {Ue)^ = Ri 2 { 9 i 2 )^ = R23(023)-^i2(0i2))) 

viii) the ( 0 f 2 , 6 'i 3 )-(^ 23 ) ^ 12 ) scheme (in which (i/e)^ = Ri 2 { 9 l 2 )R 2 ^{ 9 l^), Uu = R 23 ( 6 * 23 )^ 12 (^ 12 ))- 
For these two schemes results were given earlier in [^. We have updated the predictions 

using the best fit values of sin^ 0 i 2 , sin^ 0 i 3 and sin^ 623 , found in the most 


obtained in 30 


recent analyses of the neutrino oscillation data. 

We have not presented predictions for the BM (LG) symmetry form of Ui, in Tables]^ and 
because for the current best fit values of sin^ 0 i 2 , sin^ 023 , sin^ 0 i 3 the corresponding sum 
rules were found to give unphysical values of cos 5 (see, however, ref. p4|). 

We have found that the predictions for cos <5 of the 0f2 — ( 023 )^ 12 ) 0f3 — (023; 012 ) 

schemes for each of the symmetry forms of Up considered differ only by sign. The 0i2~(023) 012 ) 
scheme and the ( 0 f 2 , 0 f 3 ) — ( 023 ) 012 ) scheme with cu = 0 provide very similar predictions for 
cos S. 

In the schemes with three rotations in Up we consider, cos <5 is predicted to have values 
which typically differ significantly (being larger in absolute value) from the values predicted 
by the schemes with two rotations in Up discussed by us, the only exceptions being two cases 
(see Table 1^. 

We have found also that the predictions for cos <5 of the set-ups denoted as (012,023) “ 
( 023 , 012 ) ( 0 f 3 , 023 ) — ( 023 , 0 i 2 ) differ for each of the symmetry forms of Up considered both 

by sign and magnitude. If the best fit value of 623 were 7r/4, these predictions would differ 
only by sign. In the case of the (0i2,0f3) ~ (023)012) scheme, the predictions for cosd depend 
on the value chosen of the phase ui. 

We have performed next a statistical analysis of the predictions a) for cosd and Jcp 
using the latest results of the global fit analysis of neutrino oscillation data, and b) for cos 6 
using prospective sensitivities on the PMNS mixing angles. This was done by constructing 
likelihood functions in the two cases. 

For the reasons related to the dependence of cos J on w we did not present results of the 
statistical analysis for the ( 0 f 2 )0f3)- (023)0^2) scheme. This can be done in self-consistent 
models of neutrino mixing, in which the value of the phase cu is fixed by the model. 

We have found also that in the case of the 0f2 — (023, 0i2) scheme, the results for as a 
function of S or Jcp are rather similar to those obtained in in the ( 0 i 2 , 023 ) “ ( 023 ) 0 i 2 ) 
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up. The main difference between these two schemes is the predictions for sin^ 623 , which can 
deviate only by approximately 0.5 sin^ 613 from 0.5 in the first scheme, and by a significantly 
larger amount in the second. Similar conclusions hold comparing the results for the — 
(^ 23 )^ 12 ) scheme and in the “ (^ 23 )^ 12 ) scheme. Therefore, in what concerns these 

four schemes, given the above conclusions and the fact that for the {OI 2 , ^ 23 )“ (^ 23 > ^ 12 ) scheme 
detailed results already exist in the literature (see ), we have presented results of statistical 

analysis of the predictions for cos S and the Jcp factor only for the ( 6 *^ 3 , 0 I 3 )-m3^0^12) scheme. 
This was done for the five symmetry forms considered — TBM, BM (LC), GRA, GRB 
and HG. We have found, in particular, that for a given symmetry form, cos 5 is predicted 
to have opposite sign to that predicted in the { 6 I 2 , 013 )- (^23, ^^ 2 ) scheme. Thus, in the 
(^ 13 )^ 23 ) “ (^ 23 )^ 12 ) scheme analysed in the present article, one has cos<5 > 0 in the TBM, 
GRB and BM (LC) cases, and cos <5 < 0 in the cases of GRA and HG symmetry forms of U^. 
As in the (^ 12 ) ^ 23 ) “ (^ 23 )^ 12 ) set-up, there are significant overlaps between the predictions 
for cos (5 for the TBM and GRB forms, and for the GRA and HG forms, respectively. The 
BM (LC) case is disfavoured at more than 2a confidence level. Due to the fact that the best 
fit value of sin^ 023 < 0.5, the predictions for cos (5 for each symmetry form, obtained in the 
discussed two set-ups differ not only by sign but also in absolute value. We found also that in 
the ( 0 f 3 , 023 ) “ (^ 23 ’^ 12 ) scheme relatively large CP-violating effects in neutrino oscillations 
are predicted for all symmetry forms considered, the only exception being the case of the BM 
symmetry form. 

In the case of the 9^2 — (023) ^ 13 ) ^ 12 ) 0^3 — (023, 0 i 3 , 012 ) schemes we have performed 

statistical analyses of the predictions for cos <5 and the Jqp factor for the five sets of values 
of the angles [0i3,0i2] listed in Tables [3| and These sets differ for the two schemes. For 
the values of [ 0 i 3 , 0 i 2 ] given in Tables^ and H the allowed intervals of values of sin^ 0 i 2 in 
the two schemes, in particular, satisfy the requirement that they contain the best fit value 
and the 1.5a experimentally allowed range of sin^ 0i2. In the discussed two schemes the value 
of sin^ 023 is determined by the values of ^13, 0^3 and 023 (see Table . In the statistical 
analyses we have performed 623 was set to (—vr/d). Setting sin^ 0i3 to its best fit value, in the 
scheme 0f2 — (^ 23 )^ 13 )^ 12 ) ^13 ~ 0, 7r/20, vr/lO and sin“^(l/3) we found, respectively: 

sin^ 023 = 0.488, 0.501, 0.537 and 0.545. For the same values of sin^ 0i3 and 033 we obtained 
in the scheme 0f3 — ( 023 )^ 13 )^ 12 )- sin^ 023 = 0.5 1 2 , 0.499, 0.463 , 0.455. 

Further, the statistical analyses we have performed showed that for each of the two 
schemes, the five cases considered form two groups for which cosd differs in sign and in 
magnitude (Figs. and [^. This suggests that distinguishing between the two groups for each 
of the two schemes considered could be achieved with a not very demanding (in terms of 
precision) measurement of cos 5. In the analyses performed using the prospective sensitivities 
on sin^ 0 i 2 , sin^ 0 i 3 and sin^ 023 , assuming the current best ht values of sin^ 0 i 2 , sin^ 0 i 3 will 
not change, we have chosen as potential best fit values of sin^ 623 those predicted by the two 
schemes in the five cases considered (the values are listed in the preceding paragraph). These 
analyses have revealed, in particular, that for each of the two schemes, distinguishing between 
the cases inside the two groups which provide opposite sign predictions for cos <5 would be 
more challenging in terms of the requisite precision on cos <5; for certain pairs of cases pre¬ 
dicting cos (5 < —0.5 in the scheme 0f3 — ( 023 )^ 13 )^ 12 )) this seems impossible to achieve in 
practice. These conclusions are well illustrated by Figs. and However, we have found 
that, depending on the chosen potential best fit value of sin^ 023 , some of the cases are strongly 
disfavoured. Thus, a high precision measurement of sin^ 623 would certainly rule out some of 
(if not all) the cases of the two schemes we have considered. 
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The analysis performed of the predictions for the Jcp factor showed that in the 9\2 ~ 
(023)^13)^12) set-up, the CP-conserving value of Jcp = 0 is excluded at more than Su with 
respect to the confidence level of the corresponding minimum, in two cases, namely, for 
[0J'3,0 j' 2] = K/20, —7r/4], [7r/20,vr/G] (denoted in the text as Cases III and V). In the other 
three cases in spite of the relatively large predicted best fit values of Jcp, Jcp = 0 is only 
weakly disfavored (Fig. |^. For the 0f3 — (023)^13)^12) scheme, Jcp = 0 is excluded at more 
than 3cJ (with respect to the confidence level of the corresponding minimum), only in one 
case (denoted as Case I in the text), namely, for [0i3,0i2] = [ 71 / 20 , 71 / 4 ] (Fig. [^. 

The results obtained in the present article confirm the conclusion reached in earlier similar 
studies that the measurement of the Dirac phase in the PMNS mixing matrix, together with 
an improvement of the precision on the mixing angles 0 i 2 , 0 i 3 and 023 , can provide unique 
information as regards the possible existence of symmetry in the lepton sector. These mea¬ 
surements could also provide an indication about the structure of the matrix Ue originating 
from the charged lepton sector, and thus about the charged lepton mass matrix. 
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A Appendix: Statistical Details 


In order to perform a statistical analysis of the schemes considered we use as input the latest 


results on sin^ 0 i 2 , sin^ 0 i 3 , sin" 
oscillation data performed in 


023 and 6 , obtained in the global analysis of the neutrino 


15 . The aim is to derive the allowed ranges for cos <5 and JcP) 


predicted on the basis of the current data on the neutrino mixing parameters for each scheme 
considered. Eor this purpose we construct the function in the following way: = 

with Xi = {sin^ 0i2, sin^ 0i3, sin^ O 23 , 5}. The functions xf have been extracted from 
the 1-dimensional projections given in 15 and, thus, the correlations between the oscillation 


parameters have been neglected. This approximation is sufficiently precise since it allows one 
to reproduce the contours in the planes (sin^ 023 , 6 ), (sin^ 0 i 3 , 6 ) and (sin^ 023 , sin^ 0 i 3 ), given 
, with a rather high accuracy (see (3^). We construct, e.g., x^(cos J) by marginalising 


m 

2 
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over the free parameters, e.g., sim 0i3 and sim 023, for a fixed value of cosJ. Given 
the global fit results, the likelihood function. 


L(cos (5) oc exp ( — 


X^(cos(5) 


(114) 


represents the most probable values of cosJ in each considered case. When we present the 
likelihood function versus cosJ within the Gaussian approximation we use Xq ~ YliiVi — 
yi)‘^lay., with yi = {sin^ 0 i 2 , sin^ 0 i 3 , sin^ 023 }, are the potential best fit values of the 
indicated mixing parameters and ciy. are the prospective Icr uncertainties in the determination 
of these mixing parameters. More specifically, we use as Icr uncertainties i) 0.7% for sin^ 0i2, 
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which is the prospective sensitivity of the JUNO experiment 49 , ii) 5% for sin^ 023) obtained 


from the prospective uncertainty of 2% on sin^ 2023 expected to be reached in the NOvA 
and T2K experiments, and iii) 3% for sin^0i3, deduced from the error of 3% on sin^ 20i3 
planned to be reached in the Daya Bay experiment [3,51 . 


B Appendix: sin^6*23 in the (6*13,6*23) “ (^235^12) Set-up 

For completeness in Fig. we give as a function of sin^ 023 for the scheme 

The dashed lines represent the results of the global fit [15|, while the 







Figure 8: as a function of sin^ 023 in the scheme ( 0 ^ 3 , 0 I 3 The dashed 

lines represent the results of the global fit [^, while the solid ones represent the results we 
obtain for the TBM, BM (LC), GRA (upper left, central, right panels), GRB and HG (lower 
left and right panels) neutrino mixing symmetry forms. The blue (red) lines are for the NO 
( 10 ) neutrino mass spectrum. 
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solid ones represent the results we obtain for each of the considered symmetry forms of the 
matrix Ui,, minimising the value of for a fixed value of sin^ 023- The blue lines correspond 
to the NO neutrino mass spectrum, while the red ones are for the 10 one. 


References 

[1] K. Nakamura and S. T. Petcov, in K. A. Olive et al. (Particle Data Group), Chin. Phys. 
C 38 (2014) 090001. 

[2] S. K. Agarwalla et al. [LAGUNA-LBNO Collaboration], JHEP 1405 (2014) 094; 
C. Adams et al. [LBNE Collaboration], arXiv:1307.7335 [hep-ex]. 

[3] A. de Gouvea et al. [Intensity Erontier Neutrino Working Group Collaboration], 
arXiv: 1310.4340 [hep-ex]. 

[4] C. Giunti and M. Tanimoto, Phys. Rev. D 66 (2002) 113006; see also: C. Giunti and 
M. Tanimoto, Phys. Rev. D 66 (2002) 053013. 

[5] Z. z. Xing, Phys. Lett. B 533 (2002) 85. 

[ 6 ] P. H. Frampton, S. T. Petcov and W. Rodejohann, Nucl. Phys. B 687 (2004) 31. 

[7] S. T. Petcov and W. Rodejohann, Phys. Rev. D 71 (2005) 073002. 

[ 8 ] A. Romanino, Phys. Rev. D 70 (2004) 013003. 

[9] C. Duarah, A. Das and N. N. Singh, arXiv:1210.8265 [hep-ph]; S. Gollu, K. N. Deepthi 
and R. Mohanta, Mod. Phys. Lett. A 28 (2013) 31, 1350131; G. Altarelli, F. Feruglio and 
1. Masina, Nucl. Phys. B 689 (2004) 157; K. A. Hochmuth, S. T. Petcov and W. Rode¬ 
johann, Phys. Lett. B 654 (2007) 177; C. H. Albright and W. Rodejohann, Phys. Lett. 
B 665 (2008) 378. 

[10] W. Chao and Y. j. Zheng, JHEP 1302 (2013) 044. 

[11] Y. Shimizu and M. Tanimoto, Mod. Phys. Lett. A 30 (2015) 01, 1550002. 

[12] L. J. Hall and G. G. Ross, JHEP 1311 (2013) 091; Z. Liu and Y. L. Wu, Phys. Lett. B 
733 (2014) 226; S. K. Garg and S. Gupta, JHEP 1310 (2013) 128. 

[13] C. H. Albright, A. Dueck and W. Rodejohann, Fur. Phys. J. C 70 (2010) 1099. 

[14] S. M. Bilenky, J. Hosek and S. T. Petcov, Phys. Lett. B 94 (1980) 495. 

[15] F. Capozzi et al., Phys. Rev. D 89 (2014) 093018. 

[16] M. C. Gonzalez-Garcia, M. Maltoni and T. Schwetz, JHEP 1411 (2014) 052. 

[17] S. T. Petcov, Phys. Lett. B 110 (1982) 245. 

[18] F. Vissani, hep-ph/9708483; V. D. Barger, S. Pakvasa, T. J. Weiler and K. Whisnant, 
Phys. Lett. B 437 (1998) 107; A. J. Baltz, A. S. Goldhaber and M. Goldhaber, Phys. 
Rev. Lett. 81 (1998) 5730. 


35 



[19] P. F. Harrison, D. H. Perkins and W. G. Scott, Phys. Lett. B 530 (2002) 167; P. F. Har¬ 
rison and W. G. Scott, Phys. Lett. B 535 (2002) 163; X. G. He and A. Zee, Phys. Lett. 
B 560 (2003) 87; see also: L. Wolfenstein, Phys. Rev. D 18 (1978) 958. 

[20] Y. Kajiyama, M. Raidal and A. Strumia, Phys. Rev. D 76 (2007) 117301. 

[21] L. L. Everett and A. J. Stuart, Phys. Rev. D 79 (2009) 085005. 

[22] W. Rodejohann, Phys. Lett. B 671 (2009) 267; A. Adulpravitchai, A. Blum and W. Rode- 
johann. New J. Phys. 11 (2009) 063026. 

[23] J. E. Kim and M. S. Seo, JHEP 1102 (2011) 097. 

[24] J. Gehrlein, J. P. Oppermann, D. Schafer and M. Spinrath, Nucl. Phys. B 890 (2014) 
539. 

[25] A. Meroni, S. T. Petcov and M. Spinrath, Phys. Rev. D 86 (2012) 113003. 

[26] D. Marzocca, S. T. Petcov, A. Romanino and M. Spinrath, JHEP 1111 (2011) 009. 

[27] S. Antusch, C. Gross, V. Maurer and G. Sluka, Nucl. Phys. B 866 (2013) 255. 

[28] M. C. Ghen and K. T. Mahanthappa, Phys. Lett. B 681 (2009) 444; M. C. Chen, 
J. Huang, K. T. Mahanthappa and A. M. Wijangco, JHEP 1310 (2013) 112. 

[29] 1. Girardi, A. Meroni, S. T. Petcov and M. Spinrath, JHEP 1402 (2014) 050. 

[30] S. T. Petcov, Nucl. Phys. B 892 (2015) 400. 

[31] D. Marzocca, S. T. Petcov, A. Romanino and M. C. Sevilla, JHEP 1305 (2013) 073. 

[32] E. Capozzi et al, arXiv:1312.2878vl [hep-ph]. 

[33] P. 1. Krastev and S. T. Petcov, Phys. Lett. B 205 (1988) 84. 

[34] 1. Girardi, S. T. Petcov and A. V. Titov, Nucl. Phys. B 894 (2015) 733; see also: 1. Gi¬ 
rardi, S. T. Petcov and A. V. Titov, Int. J. Mod. Phys. A 30 (2015) 1530035. 

[35] S. E. King and C. Luhn, Rept. Prog. Phys. 76 (2013) 056201. 

[36] E. Plentinger, G. Seidl and W. Winter, Nucl. Phys. B 791 (2008) 60; W. Winter, Phys. 
Lett. B 659 (2008) 275. 

[37] S. Niehage and W. Winter, Phys. Rev. D 78 (2008) 013007. 

[38] S. Antusch and S. P. King, Phys. Lett. B 631 (2005) 42. 

[39] S. Antusch and V. Maurer, Phys. Rev. D 84 (2011) 117301. 

[40] S. E. King, T. Neder and A. J. Stuart, Phys. Lett. B 726 (2013) 312. 

[41] C. Hagedorn, A. Meroni and E. Molinaro, Nucl. Phys. B 891 (2015) 499. 

[42] C. Luhn, Nucl. Phys. B 875 (2013) 80. 


36 



[43] G. Altarelli et al., JHEP 1208 (2012) 021; G. Altarelli, F. Feruglio and L. Merlo, Fortsch. 
Phys. 61 (2013) 507; F. Bazzocchi and L. Merlo, Fortsch. Phys. 61 (2013) 571. 

[44] F. Bazzocchi, arXiv: 1108.2497 [hep-ph]. 

[45] W. Rodejohann and H. Zhang, Phys. Lett. B 732 (2014) 174. 

[46] R. d. A. Toorop, F. Feruglio and C. Hagedorn, Phys. Lett. B 703 (2011) 447; G. J. Ding, 
Nucl. Phys. B 862 (2012) 1. 

[47] S. F. King, C. Luhn and A. J. Stuart, Nucl. Phys. B 867 (2013) 203. 

[48] P. Ballett et al, JHEP 1412 (2014) 122. 

[49] Y. Wang, PoS Neutel 2013 (2013) 030. 

[50] P. Goloma, H. Minakata and S. J. Parke, Phys. Rev. D 90 (2014) 9, 093003. 

[51] G. Zhang [Daya Bay Collaboration], arXiv:1501.04991 [hep-ex]. 


37 



